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ABSTRACT 

Galaxy-galaxy lensing has emerged as a powerful probe of the dark matter halos 
of galaxies, but is subject to contamination if intrinsically aligned satellites of the 
lens galaxy are used as part of the source sample. We present a measurement of 
this intrinsic shear using 200,747 lens galaxies from the Sloan Digital Sky Survey 
(SDSS) spectroscopic sample and a sample of satellites selected using photometric 
redshifts. The mean intrinsic shear at transverse separations of 30-446/z -1 kpc is 
constrained to be —0.0062 < A7 < +0.0066 (99.9 per cent confidence, including 
identified systematics), which limits contamination of the galaxy-galaxy lensing signal 
to at most ~ 15 per cent on these scales. We present these limits as a function of 
transverse separation and lens luminosity. We furthermore investigate shear calibration 
biases in the SDSS, which can also affect galaxy-galaxy lensing, and conclude that the 
shear amplitude is calibrated to better than 18 per cent. This includes noise-induced 
calibration biases in the ellipticity, which are small for the sample considered here, but 
which can be more important if low signal-to-noise or poorly resolved source galaxies 
are used. 
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1 INTRODUCTION 

Weak gravitational lensing has emerged as a powerful tool 
for directl y measuring the matter distribution in the un i- 
vcrsc (e.g. iBartelmann fc Schneider! 1200 it iRefreeierl 12003). 
One of its applications has been the measurement of the 
projected matter d e nsity of g alaxies llBrainerd et al.lll996t 
Hudson et all | 1998| : iFischer et aT] 120001: ISmith et allboOli 
McKay et all I2OO it iGuzik fc Seliakl l2002t IHoekstra et alJ 
20031: ISheldon et alJ 120041: IHoekstra et alJ l2004h. groups, 
and clusters jHoeksfaa^taLf l200lt IHoekstra et alJ Il99ct 
ISheldon et al.ll200f[) . In these measurements, the observable 
quantity is the tangential shear j-t of distant source galax- 
ies induced by the lens. It is however possible that some 



of the "source" galaxies whose shears are used to measure 
lensing may in fact be physically associated with the lens; 
in this case they may be intrinsically aligned with the lens, 
producing a false lensing signal. This false signal may be 
quantified by the "intrinsic shear" A7, which is the spuri- 
ous estimate of the tangential shear that would be obtained 
by applying a shear estimator to a population of physically 
associated galaxies (see Appendix [X] for a more precise def- 
inition). This intrinsic shear is essentially a type of galaxy 
density-shear correlation, i.e. it measures the correlation be- 
tween the shears of some galaxies and the positions of oth- 
ers; it is distinct from the shear-shear correlations that con- 
taminate the gravitational shear autocorrelation measured 
by lensing surveys in the field I Van Wacrbeke^^L 2000; 
iBacon et al]l200d: iRhodes et afll200 it IHoekstra et aljhooa 
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[ Van Waerbeke et allliooS Ijarvis et alJ 120031: iBrown et all 

Several investigators have in vestigated t h is po tential 
contaminant of the lensing signal. iLee fc Pen! J200ll) mea- 
sured the density-shear correlation using 12,122 spiral galax- 
ies. These authors were primarily interested in the rotation 
axes of galaxies and their possible use as tracers of the large- 
scale tidal field (and hence the cosmic density perturba- 
tions; (Le^^^e^loO^), and so they measure position angles 
rather than shear. A crude (and model-dependent) conver- 
sion of their measurements into an intrinsic shear suggests 
A7 = -0.0037 ± 0.0025 at separations r ZD = 500/iT 1 kpc 
(see Appendix^. 

iBernstein fc Norberd l|2002t ) used 1,819 satellite galax- 
ies selected spectroscopically fro m the Two Degree Field 
(2dF) Galaxy Redshift Survey dColless et alJ 1200 ill with 
ellipticities me asured in the Autom atic Plate Measuring 
(APM) survey dMaddox et alJll99Ch to set an upper limit 
on the intrinsic shear of the satellite galaxies of IA7I ^ 0.01 
(at 95 per cent confidence) 1 for primary-satellite pairs sep- 
arated by a transverse distance of < 350/i _1 kpc. However, 
this upper limit is not presented as a function of trans- 
verse separation or lens luminosity. Circularization of the 
galaxy images due to the APM point-spread function (PSF) 
may weaken the upper lim it slightly (by ~ 10-20 per cent; 
IBernstein fc NorberdfeoO^ but this does not alter the con- 
clusion that intrinsic alignment co ntamination of the ga laxy- 
galaxy lensing signa l observed bv lFischer et alJ i200Cf) and 
iMcKav et, a?T(|200ll) is limited to < 20 per cent. 

There have been other observational studies of intrin- 
sic alignments, including reported detectio ns of higher-order 
density-shear statistics llLee fc Penll2002l) and s hear-shear 
correlations (iPen et alJ l2000t IBrown et alJl2002fl : unfortu- 
nately there is no simple conversion from these measure- 
ments into an intrinsic shear A7. 

Theoretical work on galaxy intrinsic alignments has 
been motivated both by inherent interest and (more re- 
cently) by contamination of cosmic shear surveys for which 
the intrinsic-alignment parameter of interest is the shear- 
shear correlation instead of the density-shear correlation. 
A number of estimates of the shear-shear correlations 

- both in TV-body simulations and in an a lytical models 

- have been made llCroft fc Mctzlcr 2000; Heavens et 
200d:lLee fc Penl2000tlLee fc PeJ200ll:ICatelan et alJl200ll 
Crittenden et all200ll:|Hui fc Zhandl2002l: |jindl2002lh Usu- 
ally the galaxy is approximated as a disk aligned perpen- 
dicular to the halo angular momentum vector (for spi- 
rals) or as an ellipsoid homologous with the ellipticity of 
its ha lo (for ellipticals), al though comparison to observa- 
tions (|Hej3rrans_^t L ^lJ|2i)0^) a nd sim ulations of galaxy for- 
mation iyVMi^^n^oscIi^taTI l2002t ) suggest that this pic- 
ture is too simple and overestima tes the shea r -shea r cor- 
relations. Of these authors, only ILee fc Pen! <l200ll) and 
iHui fc Zhand <l2002l) estimate the density-shear correlation 
(and even then onl y for spirals). Usi ng the rough conversion 
of Appendix^ the lLee fc Pen! <l200ll) prediction corresponds 
to A7 = —0.004 for the spiral sources , the — sign indicating 
radial alignment. IHui fc Z hang (2002) present some predic- 



1 Bernstein & Norbcrg (2002) reported an upper limit on the el- 
lipticity of 0.02, which we have converted here to a shear. 



tions for the two-dimensional (projected) density-ellipticity 
correlation as a function of angular separation, but do not 
provide a numerical estimate A7(r). 

In this paper, we use SDSS spectroscopic galaxies 
(as lenses) and photometric galaxies (as sources) in a 
luminosity-dependent study of the shear around lens galax- 
ies due to intrinsic alignments. The large number of galaxies 
in the SDSS photometric sam ple allows tighter con s traint s 
on A7 than were obtained bv IBernstein fc Norberd (2002). 
It also provides the statistical power to compute these upper 
limits as a function of transverse separation and lens lumi- 
nosity. To make the interpretation for weak lensing as simple 
as possible, we directly compute A7 using a lensing estima- 
tor (including P SF corrections) ra ther than the position- 
angle statistics of lLee fc Pen! i200 ll) or the unweighted mo- 
ments computed by the APM survey. 

Correct calibration of the shear estimator is necessary 
in order to interpret either an intrinsic alignment or a weak 
lensing-induced shear, particularly if the shear is detected 
at high signal-to-noise. Errors in shear calibration can come 
from a variety of effects such as incomplete PSF correc- 
tion, selection effects, noise-induced biases, uncertainties in 
PSF reconstruction, and incomplete knowled ge of the el- 
lipticity distribution of the source population IIKaiseill2b00l : 
IBernstein fc Jarvisi2o61lHirata fc Seliakl2003l) . Since we do 
not detect an intrinsic alignment signal, the only effect of 
calibration uncertainty is a minor degradation of our up- 
per limits on A7. Nevertheless, the calibration results pre- 
sented here are of direct interest for ongoing galaxy-galaxy 
weak lensing studies with the SDSS; calibration errors at 
the ^16-18 per cent level found here are potentially impor- 
tant in those cases where the lensing signal is detected at 
> 6(J. The noise-induced calibration biases have not been 
explicitly studied previously; while we find that this effect 
is insignificant for our sample of relatively bright (r < 21) 
galaxies, the calibration bias can be large for galaxies de- 
tected at low signal-to-noise ratio, especially if they are also 
not well- resolved. 

This paper is organized as follows: Sec. |2] includes a 
description of the lens catalog, and Sec. ^includes discussion 
of the source catalog, including details of the apparent shear 
measurement, the associated calibration uncertainties, and 
the photometric redshifts. In Sec.|lJ we describe the methods 
used to compute uncertainties on the shear. Sec. includes 
results for intrinsic and gravitational shear, and the results 
of tests for systematic errors affecting these quantities. We 
conclude in Sec. [5] with a discussion of the implications of 
these results. The noise-induced calibration bias derivation 
is presented in Appendix ICl 

A note about the cosmological model and units used in 
this paper: we compute the proper transverse separation of 
the lens and source using the physical angular diameter dis- 
tance cLa{zi) to the lens. This distance is determined assum- 
ing a flat ACDM universe with Q m = 0.27 and Q,a = 0.73, 
in accord with cosmological parameter deter minations from 
the W ilkinson Microwave Anisotropy Probe dSpergel et alJ 
2003). We present distances in kpc, where h is the re- 
duced Hubble parameter: Ho = WOh km/s/Mpc. When ap- 
propriately scaled by h, the angular diameter distance is 
only weakly sensitive even to extreme variations in cosmol- 
ogy: for a lens at redshift z — 0.3 (among the most distant 
in our sample; see Fig. [SJ, we find that cLa{zi) is reduced 
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by 7 per cent if we switch from ACDM to an open universe 
(f2 m = 0.27, £7a = 0), and by 12 per cent if we switch to an 
Einstein-de Sitter universe (Q m = 1, = 0). Given that 
we see no detection of intrinsic alignments, the dependence 
of our results on fi m and Qa is negligible for the purpose of 
constraining contamination of the weak lensing signal. 



2 LENS CATALOG 

The "lenses" used in this investigation are obtained from the 
SPSS main galaxy spectroscopic sample lYork^et_alJ l2000|; 
[ Strauss et ail 120021: iBlanton et alJ l2003at lAbazaiian et alJ 
2004). Our sample contains 200,747 galaxies within a 
solid angle of 2,499 square degrees (SDSS Sample 12; 
IBlanton et"all l2003al i. This sample is appro ximately flux- 
limit ed at Petrosian magnitude r < 17.77 iStrauss et al 
2002). (The r filter is centered at 622 nm: iFukugita et al 



Lens luminosity distribution in r 



1996). We have taken only objects at redshift z > 0.02; at 
lower redshifts, the summation of pairs of galaxies out to 
l/i -1 Mpc transverse separation becomes computationally 
expensive, and little is gained anyway because of the small 
inverse critical surface density 2 E" 1 . These spectra have 
been proce ssed by an independ ent spectroscopic pipeline at 
Princeton ijSchleeel et al.1 2004); comparisons with the sur- 
vey pipeline show that the failures occur < 1 per cent of 
the time. The o bjects have be e n ext inction-corrected us- 
ing the model of ISchleael et al.l il998h extrapolated to the 
SDSS filters using an extinction-to-reddening ratio Rv = 
3.1, and if-corrected to z = 0.1 using kcorrect vl_ll 
iBlanton et alJl2003bl) . 

Fig. shows the magnitude distribution of the lens 
galaxies, and Fig. |5] shows their redshift distribution, over- 
all and as a function of luminosity. We have used model 
magnitudes (i.e. magnitudes deter mined from the best-fit 
de Vacouleurs or exponential profile; IStouehton et al.ll2002T) 
for both the lens and the source galaxies in this paper. 



3 SOURCE CATALOG 

3.1 Apparent shear measurement 

In this section we discuss the methodology for measuring 
the apparent shear of source galaxies. In general, this must 
consist of the following steps: 

(i) Source galaxy detection and selection; 

(ii) Point-spread function (PSF) determination; 

(iii) Measurement of galaxy ellipticity, including PSF cor- 
rection; 

(iv) Conversion of ellipticity measurement into a shear. 

If we intend to measure the radial profile of galaxy halos, 
we must also use redshift and cosmological information to 
convert the shear measurements into constraints on the mass 
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Figure 1. The magnitude distribution of lens galaxies 
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Figure 2. The redshift distribution of lens galaxies, overall and 
as a function of luminosity. Note that the brighter (smaller M r ) 
subsamples peak at higher redshift. 



distribution. However, in this paper we are only interested in 
determining the spurious contribution to the apparent shear 
from intrinsic alignment of source and lens galaxies. There- 
fore we only work through the four steps outlined above. 
Also, the reconstruction of the PSF from stellar images is 
carried out by standard SDSS software dStoughton et all 
2002: lLupton et aill2004F) and will not be described here. 



2 The inverse critical surface density is denned by SJT 1 = 
<±ttGDolDls/c 2 Dos> where Dql> AlSi an d Dos are the an- 
gular diameter distances from the observer to the lens, the lens to 
the source, and the observer to the source. This quantity is useful 
because it packages the geometrical aspects of lensing together; 
the shear signal for a given weak lens at fixed physical transverse 
separation r scales 7 <x EJT . 



3.1.1 Adaptive moments 

We begin by defining the adaptive moments, which are used 
in the source selection and ellipticity measurement. The 
adaptive second moment matrix M of a galaxy's image in- 
tensity I(x, y) is found by minimizing the energy functional: 
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E(A,x ,y ,M;I) 



I(x,y) - Ae 



V/2 



dxdy, 



where: 

p 2 = (r - r ) T M _1 (r - ro) 



(1) 



(2) 



and r = (x,y). The vector ro is the object centroid, and 
the moment matrix M is taken to be symmetric and can be 
parameterized as: 



T 

M = - 



l + e + e x 
e x l-e+ 



(3) 



Here T is called the trace, and e = (e+,e x ) form the spin-2 
ellipticity. As defined here the ellipticity is restricted to lie 
in the unit circle e\ + e x ^ l. 3 The trace is one measure of 
the "size" of an object; another is the geometric mean a of 
the semi major and semi minor axes: 



2 T rv 

where e 2 = < 



Vdet M, 



(4) 



+ e 2 x . The energy functional is minimized by 



requiring: 

J(r — ro)(r — r ) T I(x, y)e~ p ^ 2 dxdy 



M 



J I{x,y)e~ p2 / 2 dxdy 



(5) 



if the centroid were also determined adaptively, we would 
simultaneously solve 



0=1 (r — ro)I(x, y)e 9 ^ 2 dxdy. 



(6) 



However the SDSS photometric pipeline (photo) fixes ro at 
the object finder's centroid rather than iteratively solving 
Eq. ©. We also use the radial fourth moment 0,4 defined 
by: 



2(1 + a 4 ) 



/ p 4 I{x,y)e- p2/2 dxdy 
J I(x,y)e~P 2 / 2 dxdy 



(7) 



Photo computes the adaptive moments M and 04 for both 
the observed (i.e. not deconvolved) galaxy image I(x, y) and 
the reconstructed PSF P(x, y). We will denote moments cor- 
responding to the galaxy image with a superscript and 
those corresponding to the PSF with a superscript ' p - ) . The 
(raw) resolution factor R2 is defined by: 



R 2 = 1 



T (P) 

TUT' 



(8) 



in the limit of a very well-resolved galaxy, R% — > 1, whereas 
for a poorly resolved galaxy the observed image is very sim- 
ilar to the PSF (/ « P) and so R 2 -> 0. 

The physical interpretation of M^ 7 ' is as a best-fit Gaus- 
sian to the image profile. The fourth moment parame- 
terizes the departure of the galaxy from Gaussianity: it is 
for a Gaussian profile, 0.17 for an exponential profile, and 
0.40 for a de Vacouleurs profile. The fourth moment of the 
PSF diagnoses the deviation of the PSF from Gaussianity. 
In particular, a PSF dominated by Kolmogorov turbulence 
(optical transfer function oc e ~ (i/i ° )5/3 : lKaiseill200fJ) should 

3 We warn the reader that the definition of ellipticity is not 
standard across the litera ture. This defi nition is consistent with 
iBernstein fc Jarvisl <2002ft . but not with lKaiser et alj <1995l) . 



have a 



(P) 



0.046, in good agreement with the typical PSF 



in the SDSS. 



3.1.2 Source galaxy selection 

Our sou rce galaxies are selected from the SDSS photometric 
catal o g jYork et al.l200(il: [ Hogg et al]l200lt [Stoughton et alj 

2002) : ISmith et al.l I2OO2I : IPier et alj booa Ubazaiian et alj 

2003) . Th e catalog is ba sed on images from the SDSS 
camera llGunn e t alj|l9 98j) processed at Princeton by the 
photo software <|Lupton et al 1l200ll:lFinkbeiner et al.ll2004l : 
iLupton et alJl20 04). In order to avoid the PSF anisot ropy- 
induced selection bias d iscussed by iKaiseJ i2000f) and 
IBernstein fc Jarvisl i2002t) . PHOTO applies a convolution to 
each image to circularize the effective PSF before running 
the object detection algorithm. We only consider objects 
that are classified by PHOTO as galaxies, are not deblended, 
do not contain saturated pixels, and do not have flags set 
indicating possible problems with the measurement of the 
image. We reject objects whose adaptive moment measure- 
ments failed, the resolution factor R$ < i, the measured 

(1)2 (1)2 

ellipticity e\ + e x > 0.95, or the radial fourth-moments 



|ai 7) | > 0.99 or 



a4 P ' I > 0.99 in at least two of the g, r, 
and i bands. Objects with extinction-corrected model mag- 
nitude fainter than 21 in r-band or with failed photometric 
redshifts (see are also rejected. In case of multiple ob- 

servations of the same object, we take the observation where 
the i-band resolution factor is largest. 

Our final source galaxy catalog contains N a =6,975,528 
galaxies, although some of these are in regions of the SDSS 
without spectroscopic coverage and hence are not actually 
used in the analysis. 



3.1.3 Measurement of galaxy ellipticities 

Since most galaxies used in weak lensing measurements 
are of comparable angular size to the PSF, the ellipticity 



(c 



CO r CO, 



determined from the galaxy image I(x,y) is gen- 



erally not equal to the true ellipticity of the galaxy. There- 
fore a correction for the PSF is necessary. The basic idea of 
any such correction is to treat the observed galaxy image / 
as a convolution of the PSF P and a pre-seeing image /: 



I{x,y)= I f(x,y')P{x~x',y~y')dxdy. 



(9) 



The problem is to determine the ellipticity of the pre-seeing 
image / from I and P. This problem has been considered by 
many author s (lKaiser_et_al . 1995; LuDDino &^£aiser[|l997; 
[ Kaiseil I2OO0I IBernstein fc Jarvisl l2002t iRefregier fc Baconl 
2003; Hirata & Scliak 2003); we use the "linear" method of 
iHirata fc S cliak (2003;) because it can be implemented using 
the PHOTO outputs. 

The "linear" method is based on the observation that if 
the galaxy and PSF are Gaussian, then the adaptive second 
moment matrix of / is related to that of the PSF and the 
pre-seeing image by: 



M 



(/) 



M (/) + M 



(p) 



implying: 



e (I) 



l)e 



(P) 



(10) 



(11) 
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Hira ta fc Seliakl ll2003l) made the approximation that the 
pre-seeing galaxy image and the PSF can be approximated 
by the quartic-Gaussian form: 



«4 , 4 



4p 2 + 2) 



-p 2 /2 



(12) 



where p is given by Eq. and similarly for P(x,y). 

Within this approximation, we can compute the corrections 
to Eq. 11011 to first order in and 04 . The calculation is 
straig htforward (although te dious) and is given in Appendix 
B of iHirata fc Seliakl (§003); it results in a correction to 
the resolution factor i?2- We compute ellipticity estimators 



(gy 



(/) zU)\ 



using the formulas presented there. Measurement 



errors on M' 1 ' are computed by PHOTO from the Fisher ma- 
trix, and propagated into errors on e^' via numerical dif- 
ferentiation of the PSF correction formulas. 

The actual image of the galaxy and the reconstructed 
PSF are both pixelized, and in principle the PSF correction 
should take this into account. Mathematically, pixelization 
is a two-step process consisting of convolution with a square 
top-hat function (i.e. integration over the pixel area), fol- 
lowed by sampling at the center of each pixel. So long as the 
convolution step is applied to both the PSF and the galaxy, 
it does not cause concern for us because Eq. @ is still valid. 
The sampling at the center of each pixel replaces the inte- 
gration in the moment computation with summation over 
pixels, which causes a problem for the integrals in Eq. 
if the integrand is insufficiently sampled (i.e. if there is a 
high-wavenumber component aliased to zero wavenumber); 
this occurs for wavenumbers: 



I > la 



2tt 



i 3.28 x 10 6 « 15.9 arcsec -1 . 



(13) 



Aliasing is more likely to occur for the PSF measurement 
than the galaxy measurement since the PSF is smaller in real 
space. If the PSF were a Gaussian with Ofwhm = 1.2 arc- 
sec (better than typical seeing in the SDSS), the function 
P(r)e~ p / 2 (where p 2 is determined from Eq.[5]using M^ p ') 
has la width of 0.36 arcsec in real space and 2.8 arcsec -1 in 
Fourier-space. Therefore, we do not expect a significant con- 
tribution to the moment integrals at the aliasing wavenum- 
ber laiias, and we have not applied any correction for pix- 
elization. 

We have determined the camera shear via diffcrcntia- 
tion of the SDSS astrometric solution <Pier et alJl2003h . and 
found it to be < 0.2 per cent in all three bands (g,r,i) used 
for the ellipticity measurement in all six camera columns. 
Given that the camera shear is small and is not correlated 
with the locations of lens galaxies, no correction for camera 
shear has been applied. Our search for a shear signal around 
random points (Sec. 15.21 produced a null result, confirming 
that camera shear contamination of the intrinsic alignment 
signal is negligible. 

The ellipticities e^' are measured separately in the g, r, 
and i bands (the u and z bands typically have lower signal- 
to- noise and so are not as useful for ellipticity measurement). 
An overall ellipticity is computed for each source galaxy by 
performing an average (weighted by the measurement error) 
of the ellipticities in the bands in which ellipticity measure- 
ment was successful. At this stage we reject objects with 
ellipticities e (/)2 > 2. 
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Figure 3. The measured RMS ellipticity of our galaxies per 
component (i.e. not the total RMS ellipticity, which is \/2 times 
greater since ellipticity has 2 components) as a function of r mag- 
nitude. The solid curve is the RMS of the measured ellipticities, 
while for the dashed curve the measurement noise has been sub- 
tracted out. For brighter galaxies, the RMS ellipticity per com- 
ponent is 0.37, which is the value we use to compute the shear 
responsivity. (The statistical significance of the downward dip at 
r ~ 20 is large, but its amplitude is only 3 per cent, which makes 
it smaller than possible calibration errors. It is therefore possible 
that this does not represent a real feature in the galaxy ellipticity- 
magnitude distribution.) 



3.1.4 Ellipticity to shear conversion 

Once the ellipticities of the source galaxies have been deter- 
mined, it is necessary to convert them into a shear estima- 
tor. The simplest such estimator is obtained by dividing the 
ellipticity by an appropriate "shear responsivity factor" 1Z: 



7 = 



TlN a 



(14) 



where N a is the number of source galaxies measured. Here 1Z 
measures the response of the mean ellip ticity to an ap plied 
shear , (e) = IZj + 0(7 2 ). It is given by: llBernstein fc Jarvisl 
2002) 



7Z — 2(1 £ r ms); 



(15) 



where e rm3 is the RMS ellipticity per component (+ or x). 
From Fig. |3] we find e rms w 0.37. (The factor of 2 comes 
from the fact that a shear of 1 per cent, when applied to a 
circular object, yields an ellipticity of 2 per cent.) 

In principle, better shear signal-to-noise ratio can be 
obtained if galaxies are wei ghted according to the ell iptic- 
ity measurement error frischer fc TvsorJl997tlMcKav et alJ 
l200ll) or if an elli p ticity-depende n t weight is used 
iHoekstra et alJ l2000l : iKaiseJ l2000l : iBernstein fc Jarvisl 
l200ab Since we only use bright (r < 21) galaxies, for which 
the measurement error is subdominant to the shape noise, 
weighting by measurement error is not helpful. Also, while 
ellipticity-dependent weighting is useful for measurement of 
gravitational shear, it complicates the interpretation of in- 
trinsic alignment studies since the apparent shear from in- 
trinsic alignments probably depends at some level on the 
method of weighting. Therefore in this paper we use only 
the "unweighted" estimator, Eq. (1141 . 

Because there is insufficient signal-to-noise ratio for a 
shear measurement around one galaxy, we report the shear 
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measurement averaged in annular bins around many lens 
galaxies. We compute the galaxy-apparent shear correlation 
function yt(r) averaged in the annular bin r_ < r < r+ by 
pair summation, 

7t(r -' r+) = A^-^)' (16) 

where the summation is over lenses a and sources j; 
N p (r~,r+) is the number of lens-source pairs with separa- 
tions between r_ and r+; and jt, a j is the + component of 
the apparent shear of source j in coordinates where the x- 
axis is aligned perpendicular to the line connecting the lens 
and source, and the y-axis is aligned along the line connect- 
ing the lens and source. That is, 7 t is the "tangential" shear 
component (so that lensing- induced shears are positive). In 
addition to 74, we report the shear 745 at a 45 degree an- 
gle to the lens (the x component in the coordinate system 
aligned with the lens-source direction) . The latter must van- 
ish in the mean by symmetry, but is nevertheless useful for 
confirmation of the uncertainty estimates. 

3.2 Calibration uncertainty 

In order to interpret the results of either intrinsic or grav- 
itational lensing shear measurements, we must understand 
two types of systematic errors in the shear measurements: 
additive biases, which are independent of the shear signal; 
and multiplicative or calibration biases, which affect the re- 
sponse of the shear estimator to an actual shear. Additive 
biases (due, e.g. to PSF ellipticity) constitute a source of 
spurious power in the shear power spectrum and hence are 
a serious issue for cosmic shear surveys that aim to mea- 
sure this power spectrum. However, since we are essentially 
computing a cross-correlation function between the galaxy 
density and shear, the additive bias can be determined by re- 
peating the galaxy-shear correlation function measurement 
using a random "lens" galaxy catalog (see Sec. 15.20 . Com- 
puting the calibration bias is much harder since there is 
no straightforward way to measure it directly. At present, 
then, the calibration bias must be understood theoretically 
by considering all known sources of calibration error. The 
major sources of error are shown in Table along with their 
estimated uncertainty. Note that these systematic error es- 
timates are in some cases very rudimentary and should be 
considered only as rough guides to the level at which sys- 
tematics might be affecting our results. 

In addition to the sources of error mentioned in the ta- 
ble, there are two other phenomena in weak lensing surveys 
that can mimic a calibration error: source redshift errors 
and stellar contamination. An error in the source redshift 
distribution results in incorrect determination of the lensing 
strength E" 1 and hence the surface density contrast AE; 
however, as our objective here is to measure apparent shear 
rather than mass, this effect is not relevant to the intrinsic 
alignment estimation. Stellar contamination of the "source 
galaxy" sample can dilute the shear signal, resulting in an 
effective calibration bias in galaxy-galaxy lensing, however 
this does not affect our measured intrinsic shear signal be- 
cause we divide our intrinsic shear signal by the fraction of 
our sources that are physically associated with the lens (as 
determined by the lens-source correlation function or "pair 
ratio;" see Sec. 01. Stars are not clustered around the lens 



Table 1. Identified sources of calibration error £7/7 in the ap- 
parent shear measurement, and rough upper limits to the level of 
error they might induce. (See text for caveats.) 



Source 


Calibration error 




(per cent) 


PSF dilution correction 


-8 to +6 


PSF reconstruction 


±1.9 


Shear selection bias 


±5.0 


Shear responsivity error 


±1.7 


Noise-rectification bias 


-1.1 to +1.4 


Total 


-18 to +16 



galaxies and so their dilution effect on the observed shear 
and the physically associated fraction cancel out in the es- 
timate of A7. 

In the remainder of this section we discuss each of the 
biases m Table □ In each case, our aim is to estimate or 
constrain the induced calibration bias £7/7, defined by: 



Sj _ 1 /9( 7+ ) | d(%) 
7 2 V d-y+ c>yx 



(17) 



T = 



where (7) is the expectation value of the shear estima- 
tor. Note that £7/7 = for an unbiased shear estimator 
((7) = 7). If ($7/7 is positive, then any galaxy-shear corre- 
lation signal will be overestimated; if 67/7 is negative, the 
signal will be underestimated. In principle, £7/7 may vary 
with angular position, however in computing the effective 
calibration bias for the galaxy-shear correlation we only care 
about the average value of Syf-y. 



3.2.1 PSF dilution correction 

As noted in Sec. 13.1.31 we do not measure the galaxy ellip- 
ticity ' directly but rather the measured image ellipticity 
e (I \ which has been diluted by the blurring effect of the 
PSF. Our correction for this effect is not perfect since the 
quartic- Gaussian form (Eq. I12H is not a perfect model for 
the ellipticity of the galaxy or PSF, and this can lead to a 
calibration error. 

The calibration error in the PSF dilution correction 
for adaptive m o ment methods was studied in detail in 
iHirata fc Seliakl j2003t) (an analysis f or the non-adaptiv e 
moment methods can be found in lErben et alJ 1200 if) . 
IHirata fc Seliakl (120031) found that for the range of param- 
eters of interest here, the "linear" correction can have cali- 
bration bias between —8 and +13 per cent (the worst case 
being exponential-profile galaxies whose resolution factor is 
near our limit R2 ~ i). However, we note that for 81 per 
cent of our sour c e gala xies, R2 *5 |, and within this range 
IHirata fc Seliakl (120031) found that £7/7 ranges from —8 to 
+4 per cent. Using this tighter constraint for the R2 ^ | 
galaxies, and the —8 to +13 per cent range for the 19 per 
cent of our sample with | ^ Ra ^ |, we constrain ^7/7 for 
the overall sample to lie in the range of —8 to +6 per cent. 
This is the range of values shown in Tabled 

It can be seen from Table Q that the PSF dilution 
correction is currently the largest item in the calibration 
error budget; this error will be reduced in future stud- 
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ies by using more accurate but computationally expen- 
sive PSF correction methods such as re-Gaussianization 
( Hira ta fc SeliakteOOSh or shao elet decomposition ijRefreigerl 



(riirata 6z acliak or shape 

l2003l:lRefregier fc Baconll2003li . 



3.2.2 PSF reconstruction 

The SDSS PSF is reconst ructed using images o f the bright 
stars by the PSP pipeline llStoughton et al1 l2002'l. Any error 
in the PSF model used in the PSF correction can translate 
into an error in the "PSF-corrected" ellipticity . Two ma- 
jor types of error concern us here: ellipticity errors and trace 
(size) errors. Ellipticity errors introduce spurious power into 
the shear measurements, which can be a serious problem 
for shear autoc orrelation measurements (IJarvis et al]l20oH: 
lHoekstrall2003) : however to lowest order these are less se- 
rious for galaxy-galaxy lensing because these errors are not 
expected to correlate with the locations of galaxies. If the 
size of the PSF is systematically over- or underesti- 

mated, this leads respectively to under- or overestimation of 
the ellipticity e^', and hence to a shear calibration error. 
We have tested the PSF traces by comparing them to the 
observed traces of "stars" (as identified by photo) with 
PSF magnitude uncertainty 4 in the range 0.05-0.1 (la). 
Fainter stars are not used because their moments are too 
noisy; brighter stars are not used because they have been 
used as part of the PSF fitting and hence the errors in their 
moments may not be representative of the typical error in 
PSF reconstruction. 

An estimate of the systematic error in the PSF trace can 
be obtained by considering the quantity q — ln(T^ 7 '/r' p '). 
If the PSF reconstruction were perfect, we would have (to 
second order in 8T W = T {1) - T (p) ): 



(9) 



(ST (1> ) 
T (p) 



1 (5T^ 2 ) 

2 T( p ) 2 



(18) 



where v is the signal-to-noise ratio, and we have used 
Eq. <C26t for the bias and variance of the adaptive trace 
estimator, assuming the ellipticity of the PSF is small com- 
pared to unity. The signal-to-noise ratio, for near-Gaussian 
PSF, is related to the PSF magnitude uncertainty apsFmag 
via v = OAlnlO/apsFmag. Thus the log trace ratio q has 
expectation value 



(q) = 0.0236 f ap SFma,A ' 



+ 0(o- PSF „ 



(19) 



which varies from 0.0059 to 0.0236 over the range of signal- 
to-noise ratios considered. 

The tail-rejected means obtained from the stars are 
(q) = -0.0087, -0.0086, and -0.0086, in g, r, and i re- 
spectively if objects more than 3a from the mean are re- 
jected; and +0.0057, +0.0057, and +0.0056 (with rejection 
at 7a). Comparison of these numbers to Eq. 11911 suggest 
that the systematic bias in q = \n(T^' /T^^) induced by 
PSF reconstruction errors is ^ 3 per cent, i.e. that if t' p ' is 
being systematically over- or under-estimated by PSP then 
the magnitude of this effect is no more than 3 per cent. 
We have therefore computed the systematic error in Tabled 



4 The PSF magnitude of an object is computed by fitting a PSF 
to the image of the object and finding the best-fit normalization. 



assuming 5T (P) /T (P) = ±0.03. The conversion to a shear 
calibration is obtained by differentiation of Eq. (111! with 
respect to T' p ', neglecting the non-Gaussian correction to 
i?2. This yields, to first order, 



8j 

7 



5e 



(/) 



e U) ~ R 2 {n2 ij T( p ) ' 



+ 



(20) 



(The error for e x is the same.) The mean value of R^ 1 — 1 
for our sample of galaxies is 0.62; multiplying by 0.03 gives 
the systematic error estimate in Table Q 



3.2.3 Shear selection bias 

If the algorithm for selecting galaxies preferentially selects 
galaxies that are nearly circular, then the estimated shear 
(as computed from Eq. 1141 will be closer to zero than the 
true shear. This could arise, e.g. because after PSF convolu- 
tion a circular galaxy is spread over fewer pixels (and hence 
will be detected at higher signal-to-noise) than a highly elon- 
gated galaxy with the same magnitude and pre-seeing area. 
The actual selection algorithm used here is more sophis- 
tocated than a simple signal-to-noise cut; we expect that 
some of our cuts, e.g. the magnitude cut at r < 21, will be 
relatively insensitive to shear selection biases whereas oth- 
ers, e.g. the resolution factor cut at R2 ^ |, will not. In 
principle the shear selection bias would be eliminated if the 
selection algorithm consisted of a simple cut on the magni- 
tude of the galaxy (which is shear-independent), however in 
practice one needs additional cuts, e.g. to remove unresolved 
galaxies from the catalog. 

The shear selection bias is difficult to compute from first 
principles because of cuts made on the galaxy sample (con- 
vergence of adaptive moments, flags, resolution factor cut, 
ellipticity cut, magnitude cut, etc.). Of these, the most wor- 
rying are the resolution factor cut and the adaptive moment 
convergence. In genera l, the shear selection bias is given by 



convergence. In genera, 
jHirata fc SeliaklboOSt 



£7+ 
7+ 



1 

n 



■InT 



(21) 



where the average is taken over galaxies used for the mea- 
surement; V is the probability of a galaxy being detected 
and selected; e is the ellipticity of the galaxy; and 7 is the 
gravitational shear. (A similar equation holds for 5*y x .) Con- 
ceptually, this equation is telling us that if the selection 
probability is decreased by applying an infinitesimal shear 
that makes the galaxy more elliptical, then the shear will be 
underestimated . 

Given the absence of a good theoretical model of selec- 
tion biases, we parameterize the problem by assuming that 
there is a selection probability that is a function of the el- 
lipticity, V(e), and that we may approximate \nV(e) using 
a polynomial of order 2a: 



c 2j e 



(22) 



3=0 



(Odd-order terms in e are forbidden by symmetry.) The 
distribution of galaxy ellipticities is shown in Fig. 21 a). In 
Fig. Elb), we have added noise to the ellipticities so that 
all the galaxies have the same ellipticity noise a e = 0.25 
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(i.e. we have added random increments Ae+ and Ae x of 
variance \/0.25 2 — a'i to each component of the ellipticity) . 
If we assume that the difference between the r < 19.0 and 
20.5 < r < 21.0 curves in Fig. Etb) is entirely due to selec- 
tion biases, then we should have: 



(a) Ellipticity distribution (measured) 



In 



n 20.5<r<21.0 (e) 



C2je 2j + const. 

3=0 



(23) 



A fit to this equation in the range < e < 1.5 using a = 4 (5 
parameters) gives C2 = —0.0178, C4 = —0.281, cq = +1.140, 
and c$ = +0.205. (We fit beyond the maximum e = 1.414 
used in this paper to suppress the characteristic "ringing" 
of polynomial fits.) Substituting this fit into Eq. 1211 . and 



noting that de/dy^ 
selection bias is: 



2e+fl 



we find that the shear 



5f 

7 



1 / f , 2 W xdlmP 
- / e( l- e) n(e)— de 



1 

7?" 



(l-e 2 Me) [J22jc 2j e 2j 



de. 



(24) 



w=i 



For the 20.5 < r < 21.0 magnitude range, we compute 
^7/7 — +0.025. We find that the computed £7/7 remains 
between 0.02 and 0.03 as we adjust the order of the fitting 
polynomial from 4 (a = 2) to 20 (a = 10), indicating that 
the polynomial fit is not affecting the results. The error given 
in Table is 5.0 per cent, which is twice the error obtained 
from Eq. 1241 (to give a more conservative estimate). We 
emphasize that we are placing a constraint on the selection 
bias, and that we have not found evidence that this bias is in 
fact present. In particular, the differences between the curves 
in Fig. Bib) could also be due to the small number of ob- 
jects with large a e (so that their ellipticity noises could not 
be "boosted" to 0.25); non-Gaussianity of the ellipticity er- 
ror distribution; systematic errors in ellipticity measurement 
(perhaps due to PSF effects); a real variation in the source 
ellipticity distribution with apparent magnitude (since the 
brighter and fainter galaxies represent distinct populations) ; 
or some combination of these effects. 



3.2.4 Shear responsivity error 

Any error in e rms translates into an error in the shear re- 
sponsivity 1Z via Eq. 1151 and hence an error in the calibra- 
tion of Eq. 1141 . While the statistical uncertainty in e rms is 
negligible because we have millions of source galaxies, the 
systematic uncertainty must be taken into account. An error 
5e rms thus translates into a calibration error given by: 



<57 

7 



STZ 
' TZ 



2 Grms $&fn 



0.86£e rr , 



(25) 



From Fig. we take Se rrns = ±0.02 as a reasonable er- 
ror estimate; if the error in e rms is larger than this then 
it must either cancel a magnitude dependence of the rms 
ellipticity or else it must affect objects of all magnitudes 
similarly (which is unlikely because the resolution factor R2 
is strongly correlated with magnitude). This leads to the er- 
ror estimate from shear responsivity error of ^7/7 = ±0.017, 
as given in Tabled 
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1.2 
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I 0.8 
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0.2 0.4 0.6 o.e 
Ellipticity e 

(b) Ellipticity distribution (with added noise, a!f=0.25) 




Figure 4. (a) The distribution of measured PSF-corrccted el- 
lipticity = \/ e \ + e x °f source galaxies as a function of 
T band magnitude. Each curve has been separately normalized 
to J n(e)de = 1. (b) The same distribution, except that each 
galaxy has had noise added to both ellipticity components so 
that the total noise is c e = 0.25 per component. This allows us 
to directly compare the various ellipticity distributions, since dif- 
ferences among them are not due to variations in <r e . A small 
number (ranging from 0.3 per cent at r < 19 to 7 per cent at 
20.5 < r < 21.0) of objects had a e > 0.25; no noise was added to 
their ellipticities. 
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Figure 5. The distribution of source galaxy resolution factors, 
dN/dR 2 = N s n(R 2 ). 
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3.2.5 Noise-rectification bias 

The fitting to determine M^'\ and the subsequent PSF cor- 
rection to obtain e' , are nonlinear operations, hence noise 
in the original image 7(r) can lead to a bias in the elliptic- 
ity e^. The additive bias (spurious power) introduced by 
noise rectification has bee n previously n oted in the context 
of co smic shear surveys iKaiseJ l200d : iBernstein fc Jarvij 
2002). Here we consider the effect of noise rectification on 
calibration bias. This effect has not been previously ana- 
lyzed, although it is present implici tly in simulations (e.g. 
iBacon et al1l200lHErben et al]l200ll) . 

The calibration correction due to noise-rectification bias 
is proportional to the noise variance and hence to v~ 2 , where 
v is the detection significance of the galaxy, i.e.: 

7 



Noise-induced calibration bias 



K N v- 2 + 0{u- 3 ). 



(26) 



The noise-rectification coefficient Kn is evaluated in Ap- 
pendix |CJ and found to be (compare to Eq. lC27l : 

K N =4(1- 3i?J 1 + R2 2 + 2e 2 rms ). (27) 

Our galaxies have resolution factor 3 < Rz < 1 and e rms = 
0.37. Within this range, Kn takes on a minimum value of 
—3.9 at R2 — I and a maximum value of +5.1 at R2 = |. A 
typical r = 21 galaxy in moderate (1.7 arcsec FWHM) seeing 
has v ~ 10, 25, and 20 in the g, r, and i bands, respectively. 
The bands are weighted by ellipticity error o^ 2 , which is 
proportional to v 2 assuming the resolution factors in the 
three bands are all equal. In this case, the weighted mean 



weighted 



Vg + Vr + V 2 ' 



(28) 



which for the signal-to-noise ratios listed above case is 
0.0027. The product Knv~ 2 appearing in Eq. 1261 is thus 
between —0.011 and +0.014. This is shown as the upper 
limit in Table □ 

While the noise-induced calibration bias (Eq. I26H is a 
subdominant source of error in this investigation, the same 
may not be true of all lensing observations. For example, 
for sources detected at 10<r (y = 10) and resolution factor 
R2 = 0.2, we find ($7/7 = +0.44, which likely dominates 
the calibration error budget for these sources (see Fig. 
The lensing studies that use these low-significance objects 
usually down-weight them due to their large ellipticity un- 
certainty <j e , so the noise-induced calibration bias in the final 
result could still be small. However in this situation the el- 
lipticity estimator e"' can be correlated with the weight, 
which leads to additional terms (positive or negative) in the 
calibration. The analysis of these correlations and their in- 
duced calibration bias is not needed here and is deferred to 
future work. 



3.3 Source photometric redshift samples 

It is critical for weak lensing experiments to have accurate 
knowledge of the source redshifts (and the lens redshifts, in 
the case of galaxy-galaxy lensing studies). It is also useful 
for studies of the intrinsic correlations, since it enables us 
to distinguish physically associated pairs from widely sep- 
arated pairs for which the gravitational lensing effect dom- 
inates the galaxy-apparent shear correlation. In this inves- 




0.4 0.6 
Resolution factor, R 2 

Figure 6. The calibration bias obtained from Eq. 1C27I for 

erms = 0.37, as a function of resolution factor R2 and detec- 
tion signal-to- noise ratio v = 5... 25. Note the large bias for poor 
resolution factors R2 < 1/3 and low signal-to- noise. The zero 
crossing at R2 ~ 0.4 is somewhat dependent on e rms and only 
exists for the first term in the asymptotic expansion of Eq. 1261 . 
hence this plot should not be taken as an indication that weak 
lensing measurements are best done with galaxies of this resolu- 
tion factor. 



tigation, we have used spectroscopic redshifts for the lens 
galaxies and photometric redshifts for the source galaxies. 
This section describes the methodology used to determine 
the photometric redshifts, and briefly considers some of tests 
of the photometric redshifts. Photometric redshift tests will 
be discussed in more detail in iMandelbaum et all J2004ft . 



3. 3. 1 Methodology 

Photometric redshift algorithms using the five SDSS band- 
passes (u, g, r, i, and z centered at 354, 4 75, 622, 763, and 
905 nm, respect ively: iFuku gita et al.lll996tl have previously 
been developed llCsabai^t^I .120031) . These photometric red- 
shift s are based on the "hybrid" photometric redshift meth- 
ods (ICsabai et al.ll200d : iBudavari et alJl2000t) . in which the 
spectral energy distribution (SED) templates are parame- 
terized; the SED parameters and source redshifts are then 
simultaneously fit to the photometric data. These are only 
available for some regions of sky (specifically the 2099 deg 2 
covered by the SDSS Data Release 1, hereafter DR1). We 
have therefore used the following methodology to derive pho- 
tometric redshifts: the DR1 galaxies in our sample are placed 
in the 5-dimensional ugriz-spa.ee. For each source galaxy 
outside of the DR1 coverage, we identify the nearest neigh- 
bour in ugriz-spa.ee with a computed photometric redshift: 



s = v/(Am) 2 + (Ag) 2 + (Ar) 2 + (Ai) 2 + (Az) 



(29) 



We report a photometric redshift failure if the nearest neigh- 
bour lies at distance s > 0.1 in ugriz-space. Of the 8,574,845 
galaxies with measured ellipticities and r < 21 that pass the 
flag cuts, matches are found for 6,925,528 (including DR1 
objects). The separation s < 0.05 for 5,889,951 of these 
galaxies, including the 3,046,566 galaxies used to construct 
the nearest-neighbour points. 
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3.3.2 Pair ratio tests 

The pair ratio R p {r) (Eq. |X2j is a useful tool for testing 
photometric redshifts; it is defined for a given transverse 
separation and lens luminosity as the ratio of the number 
density of "source" galaxies in a given redshift slice at that 
transverse separation to the number density found in the 
field. We split the lens-source pairs into widely separated 
(z 3 > zi + e, where e = 0.05 or 0.1) and nearby (|z s — zi\ < e) 
samples. If the photometric redshifts were perfect, we would 
have Rp(r) ~ 1 for the z s > zi + e samples (the exception 
to this is the magnification bias effect, which should be of 
the same order of magnitude as the shear, i.e. of order 1 per 
cent). For the \zi — z s \ < e samples, we expect R P (r) > 1 
due to galaxy clustering. 

We compute the pair ratio by dividing the number of 
lens-source pairs per unit area by the value obtained from 
randomly generated (Poisson) lens catalogs: 

n(r) 



R P (r) 



(30) 



where n(r) is the number of lens-source pairs per unit area 
in a given radial bin. The large R p for the \z\ — z s \ < e 
samples is easily seen in Fig. |7| and it is greatly suppressed 
in the z 3 > zi + e samples. The density of source galax- 
ies varies considerably with seeing because of our resolu- 
tion factor and signal-to-noise cuts. If the lens selection al- 
gorithm were perfectly independent of seeing, this would 
not contaminate the pair ratio (which is essentially a cross- 
correlation between lenses and sources), but slight varia- 
tions in the lens number density with seeing could cause 
a systematic effect (positive or negative) on Rp(r). How- 
ever, in the widely separated samples the combined effect 
of any seeing systematics and lens-source clustering (due to 
photometric redshift errors) is \R p (r) — lj < 0.03 in the in- 
nermost bins (r < 100/i -1 kpc) and < 0.01 in the outer bins 
(r ~ lOOO/i -1 kpc). When we split the lenses into different 
magnitude ranges, we find \R p (r) — 1| < 0.03 in all magni- 
tude ranges except for M r < —22, where R p ranges from 1.2 
(in 30 < r < 100ft" 1 kpc) to 1.04 (in 300 < r < 992k" 1 kpc). 
This excess is more likely to be a result of photometric red- 
shift errors than systematics in the source density because it 
is much more significant at small radii and for the brighter 
lenses. This suggests that the seeing systematic in R P {r) is 
at the 1-3 per cent level or less. 

As another test of the photometric redshifts, Fig. |H| 
shows a plot of the pair ratio for galaxies at < z s < z% — e. 
This sample should also include no physically-associated 
pairs, and therefore should have R p near 1. Instead, there 
is a large pair ratio at small radii, indicating that some 
physically-associated galaxies at higher redshift are erro- 
neously assigned to very low redshift. (We examined the pair 
ratio in the range 0.003 < z s < z\ — e to eliminate objects as- 
signed to very low photometric redshift; this did not reduce 
the pair ratio, and so we have not imposed this cut in any of 
the results presented in this paper.) Fortunately, however, 
the actual number of source galaxies represented by the ex- 
cess pair ratio in Fig. [8] is relatively small. To see this, we 
compute the excess number of galaxies in the < z 3 < zi — e, 
zi — e < z s < zi + e, and z 3 > zi + e samples: 

N+ = (TV - N rand )[z 3 > zi+e] 

N = (N - Nrand)[Zl - e < z s < zi + e] 
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Figure 7. Pair ratio for several samples, as indicated on the plot. 



7V_ 



(TV - N rand )[z s < zi-e], 



and compute the normalized fractions: 
N+.o,- 



TV+ + TV + TV_ 



(31) 



(32) 



The fractions /+ , fo, and /_ respectively represent the frac- 
tion of physically associated source galaxies whose photo- 
metric redshifts are overestimated, correctly estimated (to 
within ±e), and underestimated. For the sample of all galax- 
ies at radii 30 < r < 446/i _1 kpc, we find /+ = 0.13, 
fo = 0.82, and /_ = 0.05 (for e = 0.1), i.e. only 5 per cent 
of the source galaxies that are clustered with the lens have 
redshifts underestimated by more than 0.1. The overestima- 
tion and underestimation failure rates increase to /+ = 0.26 
and /_ = 0.16 if we narrow our redshift slice to e = 0.05. 
We find that /_ is an increasing function of lens luminosity, 
which is expected since (i) the photometric redshift can be 
an underestimate of the true redshift only if the true redshift 
exceeds e, and (ii) the more luminous lenses tend to be more 
distant. Another trend that we find is that the "probability 
of correct source photometric redshifts" fo decreases with 
increasing transverse separation r (e.g. for the all lenses, 
e = 0.1 sample we have fo = 0.90 at 30 < r < 100/i~ kpc, 
decreasing to fo = 0.80 at 300 < r < 446ft" 1 kpc), which 
could be indicative of a variation of satellite galaxy colors as 
a function of radius. The result that /+ > /_ indicates a bias 
or skewness in the photometric redshift error distribution. 

It thus appears that for the sources used in this pa- 
per, roughly 40 per cent of the source galaxies physically 
associated with a lens are scattered outside of zi ±0.05 by 
photometric redshifts errors, and 20 per cent are scattered 
outside of 2; ± 0.1. 

3.4 Excluded lens-source pairs 

At sufficiently small lens-source separation angles, the mea- 
surement of sources becomes complicated because the source 
can overlap with the lens (or with the tail of the PSF con- 
volved with the lens). Problems that occur in this regime 
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Figure 8. Pair ratio for < z s < z; — e. Counterintuitively, 
the e = 0.1 curve is above the e = 0.05 curve; this is because 
higher-redshift (and hence more luminous due to the apparent 
magnitude cut for SDSS spectroscopy) lenses are weighted more 
heavily when we examine only pairs with the "source" at lower 
redshift. 



can include incorrect ellipticity measurement if overlapping 
objects are used in the analysis, and selection bias if they are 
not. In our case deblended objects are not used for ellipticity 
measurement. 

A crude way to determine which annular bins are af- 
fected by the overlapping source-lens effect is to examine 
the pair ratio for source samples selected to be within red- 
shift \z s - zi\ < e of the lens (Fig. for e = 0.1). The 
effect of clustering is easily seen: the pair ratio is greater 
than unity and (over most of the range of radii) increases 
as the lens magnitude M r increases and the transverse sep- 
aration r decreases. However, at the innermost radii there 
is a drop off in R P (r) since sources inside these radii are 
deblended. Brighter galaxies have a more pronounced dip 
in the excess at the very innermost radii (out to 40 physi- 
cal kpc for the very brightest bin), for two reasons: (i) 
these bright galaxies are typically quite large; and (ii) they 
are more distant so the same angular separation (relevant to 
deblend ing) corresponds to g reater physical transverse sep- 
aration. |PadmaiwMwii_et_aL| i2004T) find that for M r = —22 
elliptical galaxies (the fainter edge of this brightest luminos- 
ity bin), the typical 50 per cent light radius is 10/i -1 kpc, in 
which case there may be significant overlap at the innermost 
bin (r = 20h~ 1 kpc). Since the angular diameter distance 
is 15zi kpc per arc second, and the typical seeing in the 
SDSS is 1-2 arc seconds, the deblending effect can increase 
the minimum usable radius by tens of h~ kpc for our bright- 
est subsample (typical zi ~ 0.25). For the faintest subsample 
(typical zi ~ 0.03), this increase is much less important. 

Therefore, for the brightest two luminosity samples 
(M r < —21), we have rejected the inner bins where R p is 
increasing with r. (We have verified that rejecting one extra 
bin or one fewer bin does not significantly affect the final 
results.) 



Figure 9. Pair ratio for \z s — z\\ < 0.1, several luminosity sub- 
samples. 



4 SHEAR UNCERTAINTY ESTIMATION 

We use two methods to compute the statistical uncertainty 
in the shear jt- The first is a random lens method, in which 
we use our source catalog to compute the covariance of the 
shear signal around random points within the lens sample 
region. The second is a bootstrap method in which we cut 
the lens sample region into 150 subregions that are then re- 
sampled. Because both methods use the actual source cat- 
alog, they will correctly take into account error from any 
spurious shear power in the source catalog. The random lens 
method has the disadvantage of not taking into account fluc- 
tuations in the lens number density (particularly if they cor- 
relate with the noise or systematics in the source catalog). 
The bootstrap method addresses these shortcomings on the 
smallest scales, but suffers due to the fact that the shear 
signals computed in the 150 subregions are not truly inde- 
pendent, because a source galaxy near a subregion boundary 
may contribute to the measurement of shear around lenses in 
multiple subregions. In addition, the noise in the bootstrap 
covariance matrix (which can have a significant impact on, 
e.g., x 2 testing) is more difficult to assess. The bootstrap co- 
variance matrix can be made essentially noiseless by taking 
thousands of subregions, but then the correlations among 
the subregions become large and the validity of the boot- 
strap method is undermined. 

Faced with the deficiencies of either the random lens 
method or the bootstrap method individually, we have cho- 
sen to combine them as follows. The errors on the intrin- 
sic alignment determination are computed using both the 
random-lens and bootstrap covariance matrices; consistency 
between the two methods suggests that both are reasonable 
estimates of the uncertainty. The x 2 analysis is performed 
using the covariance matrix from the random lens catalogs; 
the effect of noise in the covariance matrix is taken into ac- 
count as described in Appendix[D] The outermost bins in the 
faintest luminosity subsamples are excluded from the analy- 
sis since galaxy clustering renders the random-lens method 
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inaccurate here, and it is for these bins that the indepen- 
dence of subregions assumed by bootstrap has the least ro- 
bust justification. 

4.1 Random lens method 

We compute our statistical errors using a random lens 
test, which involves computing the shear around randomly- 
located "lenses," where Poissonian (independent) random 
lens locations were generated using the angular mask of the 
survey. We constructed 78 random lens catalogs; the galaxies 
in each random catalog were assigned magnitudes and red- 
shifts drawn from the true sample without replacement. The 
covariance matrix between radial bins i and j is computed 
as: 



Galaxy power spectrum contributions 



M 



M 

£ 



Jt(ri)yt(rj 



(33) 



where M — 78 is the number of random catalogs, the index 
a denotes the random catalogs, and 7t(?"i) is the sample 
mean of the shear in radial bin i averaged over all of the 
random catalogs. So long as the error distribution of the 
{lt(ri)} is jointly Gaussian, C gives an unbiased estimate of 
the covariance (but it cannot be used naively in \ 2 tests; see 
Appendix|HJ . The number of random catalogs used must be 
N (the number of radial bins) in order to produce a non- 
degenerate estimator C for the covariance, but is limited by 
available processor time. 

The Poisson-distributed random lens method produces 
valid error bars on angular scales where the angular (two- 
dimensional) galaxy power spectrum is Poisson-dominated 
and the clustering contribution can be safely neglected. Ne- 
glecting edge effects in the survey, this range of angular 
scales can be determined from the ratio of the galaxy clus- 
tering power spectrum to the Poisson spectrum: 



i>{k) 



P(k, clustering) 
P(k, Poisson) 



Ngal 



dN v 



dr 



-J (kr)dr, 



(34) 



where N ga i is the number of galaxies in the sample, k is 
the physical (not comoving) wavenumber of interest, r is 
the transverse separation, and N P:excess is the number of 
pairs of galaxies with transverse separation ^ r minus the 
number of pairs derived from a random (Poisson) catalog. 
We use the average of the angular diameter distances to 
the two galaxies for the conversion of angular separation 
into transverse separation r (this should be valid because no 
correlations should exist between galaxies at very different 
angular diameter distances). The ratio ip(k) is plotted for 
our samples in Fig. 1101 

If we were cross-correlating all of the lens galaxies 
against the same sample of "source" galaxies, the error vari- 
ance of the cross-spectrum would simply be re-scaled in pro- 
portion to the galaxy power spectrum, i.e. as oc l + tjj(k). (If 
survey boundary effects can be neglected it does not mat- 
ter that we are computing the correlation function via pair 
summation as opposed to via a cross-spectrum in Fourier 
space.) If ip(k) <C 1, then the Poisson power spectrum is a 
good approximation to the lens galaxy power spectrum and 
hence to the error of ft- If i>(k) 3> 1, then the lens galaxy 
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Figure 10. The ratio ip(k) of clustering to Poisson contributions 
to the galaxy power spectrum. This is obtained from the data 
using Eq. 1341 . imposing a smooth cutoff in the integrand in the 
range 20-30 h~ l Mpc. 



power spectrum is much greater than the Poisson estimate, 
and the error on j t for Poisson-distributed lenses is less than 
the error for "real" clustered galaxies. Thus we expect the 
Poisson random lens method to be a good way to estimate 
error bars in the case where ip(k) <C 1. In this paper we are 
cross-correlating each lens galaxy with the shear of source 
galaxies within the redshift slice zi ± e; this increases the im- 
portance of the clustering contribution. However, for each of 
our luminosity subsamples, the width of the lens redshift dis- 
tribution Az is of order 2e (or less, for the lowest-luminosity 
subsamples), hence ip(k) <gC 1 is still a valid condition for the 
clustering to be unimportant. Since the galaxy power spec- 
trum scales roughly as P(k, clustering) oc fc -1,2 , whereas 
the Poisson contribution is independent of k, the clustering 
contribution dominates on large scales whereas the Poisson 
contribution dominates on small scales. The random lens 
method is thus valid only at the smaller scales. 

The above reasoning applies in Fourier space. The real- 
space galaxy-apparent shear correlation function 7t(r) is re- 
lated to the galaxy-apparent shear cross-power spectrum by: 



7*M = ^ 



t (k)J 2 (kr)dlnk. 



(35) 



The J2 function peaks at kr m 3.055; therefore, we take 
only the annuli out to a distance r max = 3.055/fc m i n where 
tp(kmin) = 0.2. The values of r ma x are shown in Table|2| 

While the clustering of lens galaxies has the effect of 
increasing our error bars, the clustering of source galaxies 
around lens galaxies tends to decrease errors because there 
are more "sources" near the lenses in the real catalog than in 
the Poisson catalogs. This effect is significant when 7? p — 1 is 
of order unity (or larger), i.e. in the innermost 200 h~ kpc of 
our brighter subsamples. It is only important for the source 
samples whose redshifts are similar to the lens redshifts (the 
zi — e < z s < z\ + e samples) since there is negligible lens- 
source correlation between different redshifts. We have not 
applied any source clustering correction in the computation 
of our error bars, thus our error estimates should be viewed 
as conservative. 
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Table 2. The number of lens galaxies in each luminosity sub- 
sample, and the maximum radius r max at which the Poisson 
random-catalog uncertainties are valid. We estimate r max from 
r maI = 3.055/fc m ; n , where fc m i n is the minimum wavenumber 
(largest scale) at which ?p(k) $J 0.2, i.e. at which the clustering 
power spectrum is ^ 20 per cent of the Poisson power spectrum. 



(a) Mean shear, z. >z, +0.1 



Subsample 



gal 



fmax (h 1 Mpc) 



all 

M r < -22 
-22 < M r < -21 
-21 < Mr < -20 
-20 < M r < -19 
-19 < Mr SC -18 

M r > -18 



200747 
11337 
54768 
72464 
41276 
14475 
6427 



0.43 
>1.00 
>1.00 
0.74 
0.55 
0.52 
0.43 



4.2 Bootstrap method 

The bootstrap method, as applied here, consists of the fol- 
lowing major steps: 

(i) Split the lens catalog into M = 150 subcatalogs con- 
taining similar numbers of lens galaxies. 

(ii) Generate K — 10 5 synthesized lens catalogs. Each 
synthesized lens catalog is a concatenation of M = 150 of 
the subcatalogs, selected randomly with replacement. (Some 
lens galaxies appear in the synthesized catalog multiple 
times, others not at all.) 

(iii) Compute the shear signals {7] }«=i from each of the 
synthesized catalogs, and compute their covariance Cboot- 

In order to produce subcatalogs that are as independent 
as possible, we divide galaxies into subcatalogs that consist 
of different regions of sky. The subcatalog boundaries are 
based on the SDSS A, 77 coordinate system dStouehton et alJ 
2002). The lens galaxies are first split into 45 rings of width 
AA = 4 degrees bounded by parallels of constant A. Within 
each ring, they are then rank-ordered by rj (in the range of 
— 180 to +180 degrees). Having arranged the galaxies in this 
order, we cut the catalog into 150 approximately equal-sized 
consecutive groups, which become our subcatalogs. This re- 
sults in galaxies falling into the same subcatalog as most 
of their neighbours. (Some subcatalogs consist of multiple 
"islands" in different parts of the survey; we have not at- 
tempted any detailed optimization of the subcatalog cuts.) 

The covariance matrix computed via the bootstrap 
method is noisy even as K — > 00 because a finite num- 
ber M of subcatalogs has been used. In the extreme case 
of M ^ N, the shear signals from each of the subcatalogs 
would form an M -dimensional plane in the ^-dimensional 
shear space; then the signals computed from bootstrap re- 
sampling of these subcatalogs would all lie in this plane 
and their covariance matrix would be degenerate. As M in- 
creases, the noise in the bootstrap covariance matrix goes 
down, but the samples become less and less independent. 
We have used M = 150 here. Since some of the subcata- 
logs receive higher weight than others (due to survey ge- 
ometry, variations in source density, etc.), the calculation 
of Appendix E| cannot be directly applied to the bootstrap 
method. While this makes % 2 testing using the bootstrap 
method difficult, the errors on the final bin-averaged intrin- 




100 



1000 



r (IT 1 kpc) 

(b) Tangential shear comparison, all magnitudes 




100 1000 
r (h" 1 kpc) 

Figure 11. (a) Shear signal for z s > z\ +0.1. (b) Comparison of 
tangential shear for pairs widely separated vs. nearby in redshift. 
In both panels we have combined groups of three radial bins for 
clarity. 



sic alignments obtained from Cboot are stable and provide a 
good check of the correctness of the random-lens errors. 



4.3 Chi-squared test 

We may test the random-catalog errors using a \ 2 test, since 
the 745 shear must (in the mean) vanish by symmetry. At 
this point, it is essential to note that error covariance ma- 
trices obtained from simulations (e.g., as used in this paper) 
cannot be blindly used in place of the true covariance ma- 
trices in statistical tests, because they are obtained from 
a finite number of simulations and hence have a statistical 
error distribution themselves. In particular, when testing a 
signal for consistency with zero, it is common to use the x 2 
test, based on the variable: 



X(45) = 745 C 



745 



(36) 



Here 7 45 = (745(^1), 74s(rjv)) is the N- dimensional vec- 
tor of shears in the N radial bins. If 7 45 is Gaussian- 
distributed with covariance C and mean 0, then Xus) nas a 
X 2 distribution with N degrees of freedom. If we replace C 
with the estimate C, the distribution of x 2 becomes much 
more complicated. This distribution, while analytically in- 
tractable, is easy to compute via Monte Carlo simulations 
(see Appendix [HJ, and we have used this approach for all 
the p-values quoted in this paper. 

We test for zero 745 by computing the statistic X(4B) — 
7(45) C~ a 1 nd 7(45) , where C ra nd is the error covariance matrix 
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Table 3. The \ 2 values for fit to zero shear (for both the tangen- 
tial shear jt and the rotated shear 745) for various lens samples. 
The p values quoted take into account the effect of noise on the \ 2 
calculation as explained in Appendix [H] and represent the prob- 
ability of a random vector having \ 2 less than the given value 
by chance. The number of usable bins (i.e., radial bins where r 
is large enough to avoid overlap with the lens galaxy but small 
enough that the clustering contribution is unimportant) is N. 
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all 
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-19 


27 
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0.99998 
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all 
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26 
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0.081 


-20 < Mr < - 


19 


30 


28.3 


0.044 


21.9 


0.007 


-21 < M r < - 


-20 


32 


85.4 


0.901 


43.9 


0.225 


-22 < M r < - 


-21 


33 
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0.142 


32.6 
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27 
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0.99990 
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> 0.99999 
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0.157 








-zi\< 0.05 






all 




26 


37.0 


0.411 


25.1 


0.085 


Mr > -19 




26 


34.6 


0.331 


14.2 


0.002 


-21 < M r < - 


19 


30 


35.6 


0.149 


28.1 


0.043 


M r < -21 




33 


40.9 


0.128 


55.6 


0.424 



derived from random catalogs. (We put a "hat" on the x 2 as 
a reminder that the covariance matrix is noisy, so that this 
statistic does not have the usual x 2 distribution; see Ap- 
pendix |D|) The shear data are plotted in Fig. Illf a). while 
the X(45) values (broken down by lens luminosity) are shown 
in Table |21 along with their p- values (the cumulative proba- 
bility of obtaining a lower x 2 assuming that there is indeed 
zero signal, the errors are Gaussian, and the random real- 
izations correctly sample the error distribution). Out of the 
22 subsamples shown in the table, only 1 has a p-value for 
the 45-degree shear above 0.95, and 4 have p values above 
0.9. These numbers are consistent with the high p-values oc- 
curring due to chance. However, the high incidence of low p- 
values in the \z 3 — zi\ < e (two subsamples havep( 45 ) < 0.01) 
samples requires explanation. This is likely due to conserva- 
tive assumptions made in computing the errors, particularly 
the neglect of source-lens clustering, which leads to overesti- 
mated errors and therefore lower Xj4S) an d P~ values. Because 
source- lens clustering is more significant for the \z s — zi\ < e 
samples than for the z 3 > z\ + e samples, this explanation 
is consistent with the fact that the abnormally low p- values 



occur for the \z s — zi\ < e samples. Thus, we conclude that 
the p- values for the 745 shear show no evidence for system- 
atic effects and suggests that the random-catalog procedure 
is not detectably underestimating the error bars. 

For completeness, we also show in Table [3] the tangen- 
tial shear xft) an d P~ values. In this case, however, the signal 
is inconsistent with zero for many of the subsamples, with 
high p for the brighter luminosity lenses and widely sepa- 
rated (z s > Zi + e) lens-source pairs. As can be seen from 
the luminosity-averaged shear signal (Fig. 1111 . this is be- 
cause the measured 74 > but only for widely separated 
pairs, consistent with gravitational lensing. Since in this pa- 
per we do not fit a model to the lensing shear signal, the 
tangential xft) values are not directly useful for establishing 
the accuracy of error bars. 



5 RESULTS 

5.1 Intrinsic shear 

Table|l]shows the limits on intrinsic shear that can be com- 
puted for the various luminosity subsamples and sets of 
physical radii. The shear was computed by averaging the 
shear over the range of radii listed, weighting by the total 
number of lens-source pairs in each radial bin, 



Aj = 



N 



N-N, 



— y^w»7t(rj) = Rp y^Wi7t(n),(37) 

and '— J Hp — I 



where N and N ra nd are summed over all relevant radii 
to account for the dilution by non-physically associated 
pairs. The factor of N/(N - N ran d) = Rp/(Rp - 1) is 
used because intrinsic shear only applies to the excess pairs 
N(n) — N ra nd(ri) (see Appendix 0. This number depends 
on the set of radii and luminosities considered, and as can 
be seen from Fig. |5J is large for fainter foregrounds and for 
larger values of radii. (It ranges from close to 50 for the 
faintest subsamples, down to 1.8 for the brightest subsam- 
ples and smallest radii.) Note that shears and errors in all 
figures have not been corrected to account for this dilution 
effect, and it has no effect on the x 2 values in table[j|]because 
the dilution factor cancels out. 

The errors on A7 were computed using both random 
catalogs and bootstrap (see Sec.^J. First, the averaged shear 
was computed for each of 78 random catalog outputs (using 
the same weighting), and a' ran<i ' is the standard deviation of 
these averaged shears. This random catalog error was then 
used to compute the confidence intervals using the Student's 
^distribution, which is quite close to Gaussian for 77 degrees 
of freedom. Second, the error from the bootstrap covariance 
matrices was estimated as 



(boot) 



WiWjCboot 



(38) 



and confidence intervals were computed using a Gaussian 
distribution (since the noise in the bootstrap covariance ma- 
trix Cboot is unknown). In computing the confidence intervals 
in Table 21 we used the greater of the two errors (bootstrap 
or random catalog). 
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Table 4. Limits on the tangential intrinsic shear A7, for various subsamplcs and ranges of radii (statistical errors only). The 10 3 A7 
column shows the central (most favored) value. Quantities have been corrected to account for dilution by non-physically associated 
galaxies in the \z 3 — z t \ < e samples as described in the text. The statistical uncertainties (lcr) on A7 are shown in the next two columns 
for the random-catalog and bootstrap methods. A7 g ; is the estimated value of shear for these samples assuming a simple model for 
contamination by pairs that are actually gravitationally lensed (the errors are lcr statistical errors on shear measurement only, i.e. they 
do not include uncertainties associated with the source redshift distribution). The next two columns show the pair ratios R p for the 
\z s — zi\ < e ("in") and z s > z\ ± e ("out") samples. Finally, we show confidence intervals for A7. The "95% stat." confidence interval 
is computed using the central value of A7 — A7 9 ; and statistical errors only. The "99.9% stat.+sys." includes identified systematics as 
discussed in the text. In the row containing the [— 00, +00] confidence interval, we find no detection of physically associated lens-source 
pairs when the allowed range of systematic errors is taken into account. 
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min 


max 


min 


max 


e = 0.1 


all 




[30, 446] 


1.6 


1.0 


1.0 


1.45 ±0.17 


1.155 


1.008 


-1.8 


2.1 


-6.2 


6.6 


all 




[30, 100] 


2.7 


2.1 


1.7 


1.96 ± 0.32 


1.355 


1.007 


-3.4 


5.0 


-10.5 


12.5 


all 




[100, 300] 


-0.4 


1.1 


1.3 


1.26 ±0.22 


1.181 


1.010 


-4.1 


0.8 


-9.5 


5.1 


all 




[300, 446] 


3.6 


1.7 


1.6 


1.58 ±0.25 


1.120 


1.007 


-1.4 


5.4 


-7.4 


12.8 


M r > - 


-18 


[30, 446] 


-19.6 


14.9 


17.7 


-0.19 ±3.05 


1.022 


0.991 


-54.2 


15.3 


-612.4 


305.5 


Mr > - 


-18 


[30, 100] 


-18.1 


25.4 


22.9 


9.16 ±4.20 


1.068 


0.990 


-77.7 


23.3 


-206.6 


115.3 


Mr > - 


-18 


[100, 300] 


-24.0 


24.6 


21.4 


-3.13 ±4.34 


1.024 


0.991 


-69.8 


28.1 


-412.1 


243.1 


Mr > - 


-18 


[300, 446] 


-15.5 


24.0 


31.4 


-0.25 ±5.02 


1.017 


0.990 
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46.3 


— 00 
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-19 to 
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-0.7 
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1.013 
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— 28.6 


20.8 


-19 to 
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10.0 
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1.010 


-37.3 


2.1 
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21.5 


-19 to 
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-6.9 


5.4 
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1.8 
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1.53 ±0.33 
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-20 to 
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5.2 
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-0.7 
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1.28 ±0.49 
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1.007 


-7.1 
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1.6 


2.4 


2.3 


1.67 ± 0.45 


1.102 


1.004 


-4.8 


4.7 


-12.9 


12.8 


-21 to 


-20 


[27, 735] 


2.6 


1.3 


1.3 


1.18 ±0.20 


1.139 


1.011 


-1.1 


3.9 


-5.8 


9.8 


-21 to 


-20 


[27, 100] 


6.6 


3.5 


2.6 


2.59 ±0.57 


1.427 


1.007 


-3.0 


11.0 


-13.2 


23.2 


-21 to 


-20 


[100, 300] 


-1.0 


2.0 


1.8 


1.32 ±0.28 


1.217 


1.014 


-6.2 


1.6 


-14.1 


7.8 


-21 to 


-20 


[300, 735] 


3.5 


1.5 


1.6 


1.08 ±0.26 


1.119 


1.010 


-0.8 


5.6 


-5.8 


12.7 


-22 to 


-21 


[33, 992] 


2.8 


1.1 


1.2 


1.57 ± 0.21 


1.146 


1.006 


-1.1 


3.6 


-5.7 


8.8 


-22 to 


-21 


[33, 100] 


3.6 


3.3 


2.3 


1.80 ±0.54 


1.807 


1.027 


-4.8 


8.5 


-14.4 


19.0 


-22 to 


-21 


[100, 300] 


3.3 


1.8 


1.8 


2.00 ±0.37 


1.355 


1.019 


-2.2 


4.8 


-8.9 


12.4 


-22 to 


-21 


[300, 992] 


2.6 


1.4 


1.6 


1.48 ±0.26 


1.121 


1.005 


-2.0 


4.1 


-7.2 


10.1 


M r < - 


-22 


[45, 992] 


5.3 


1.6 


1.5 


2.25 ±0.37 


1.309 


1.048 


-0.2 


6.4 


-7.8 


18.0 


M r < - 


-22 


[45, 100] 


-0.6 


6.7 


3.9 


4.59 ± 1.28 


2.818 


1.204 


-18.4 


8.0 


-43.7 


29.8 


M r < - 


-22 


[100, 300] 


3.5 


2.8 


2.2 


2.29 ±0.53 


1.855 


1.101 


-4.4 


6.9 


-15.3 


18.8 


M r < - 


-22 


[300, 992] 


6.3 


1.9 


1.8 


2.28 ± 0.44 


1.246 


1.042 


0.2 


7.8 


-8.2 


22.1 


e = 0.05 


all 




[30, 446] 


-0.1 


0.9 


1.2 


0.61 ± 0.08 


1.201 


1.014 


-2.9 


1.6 


-7.1 


5.1 


all 




[30, 100] 


-0.7 


2.4 


1.8 


0.97 ± 0.15 


1.458 


1.022 


-6.4 


3.0 


-14.1 


9.7 


all 




[100, 300] 


-1.3 


1.2 


1.4 


0.50 ±0.10 


1.234 


1.017 


-4.6 


0.9 


-9.6 


4.5 


all 




[300, 446] 


1.5 


1.6 


1.9 


0.68 ±0.11 


1.155 


1.012 


-2.9 


4.4 


-8.5 


10.7 


M r > - 


-19 


[30, 446] 


-4.5 


4.1 


4.8 


0.51 ± 0.32 


1.067 


1.004 


-14.5 


4.4 


-29.9 


15.8 


Mr > - 


-19 


[30, 100] 


-12.2 


10.6 


9.0 


0.87 ± 0.65 


1.135 


1.005 


-34.2 


8.0 


-65.4 


31.2 


M r > - 


-19 


[100,300] 


-8.4 


5.5 


5.7 


0.00 ± 0.45 


1.076 


1.006 


-19.6 


2.8 


-38.5 


14.7 


M r > - 


-19 


[300, 446] 


1.0 


6.8 


7.8 


0.96 ±0.54 


1.055 


1.004 


-15.2 


15.3 


-36.9 


36.9 


-21 to 


-19 


[25, 545] 


-0.4 


1.2 


1.1 


0.58 ±0.08 


1.196 


1.015 


-3.3 


1.4 


-7.7 


5.1 


-21 to 


-19 


[25, 100] 


-0.8 


3.0 


2.6 


1.14 ±0.21 


1.469 


1.018 


-7.8 


4.0 


-17.0 


12.2 


-21 to 


-19 


[100, 300] 


-2.4 


1.5 


1.6 


0.50 ±0.12 


1.249 


1.018 


-6.0 


0.3 


-12.0 


4.0 


-21 to 


-19 


[300, 545] 


0.9 


1.7 


1.7 


0.57 ± 0.11 


1.163 


1.013 


-3.1 


3.7 


-8.7 


9.6 


M r < - 


-21 


[33, 992] 


2.3 


1.0 


1.0 


0.61 ± 0.07 


1.253 


1.018 


-0.4 


3.7 


-3.4 


8.0 


M r < - 


-21 


[33, 100] 


2.0 


3.8 


2.1 


1.00 ±0.24 


2.365 


1.100 


-6.5 


8.5 


-16.8 


19.5 


M r < - 


-21 


[100, 300] 


2.8 


1.7 


1.4 


0.85 ±0.13 


1.614 


1.049 


-1.4 


5.3 


-6.3 


11.6 


Mr < - 


-21 


[300, 992] 


2.1 


1.3 


1.2 


0.58 ±0.08 


1.209 


1.015 


-1.0 


4.1 


-4.7 


9.0 
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(a) Tangential shear comparison, M r >-18 



(b) Tangential shear comparison, -19<M r <-18 
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(e) Tangential shear comparison, -22<M r <-21 
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(d) Tangential shear comparison, -21<M r <-20 
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(f) Tangential shear comparison, M r <-22 
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Figure 12. Shear as a function of radial pair separation; panels (a)-(f) show the six luminosity bins in order of increasing luminosity. 
The solid line shows the widely-separated lens-source pairs (z s — z l > 0.1) and the dashed line shows the nearby pairs (|z s — Z;| < 0.1). 
Each point on the plot shows three of our radial bins averaged together to increase the signal/noise ratio for visual inspection. Note that 
the vertical scales are increased for panels (e) and (f). 



As a test of consistency, we compared the errors er (ran ' 
and a^ 00 *' computed from the random catalogs and from 
the bootstrap covariance matrices. For most of the samples 
shown in the table, the difference between the two error 
values was less than 25 per cent; the exceptions are the 
r < 100/i -1 kpc bins for the brighter lens samples, for which 
the bootstrap errors are smaller than the random catalog 
errors. This is probably due to the fact that the number 
density of source galaxies in the \z s — z\\ < e samples is 
greater near lens galaxies than in the field, hence the shape 
noise is suppressed. This effect is taken into account by the 



bootstrap but not by the random catalogs, hence the boot- 
strap errors are smaller. 

A possible source of contamination to the intrinsic align- 
ment signal is gravitational lensing of source galaxies at 
z s > z ; . We have attempted to estimate the magnitude of 
such contamination (10 3 A7 9i in the table) according to the 
following procedure. First, the mean value of the inverse 
critical surface density E~ 1 (z s , zi) for each sample was com- 
puted, where E" 1 = for z s < z\. Then, the averaged shear 
value 7 for the radii of interest was computed. The estimated 
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contamination is then 



&J 9 i = 



N 



> zi + e] \ N - N rand 



j[z s > zi + e] 



(39) 



where N/(N — N ran d) is the dilution factor mentioned pre- 
viously for the \z s — Zi\ < e samples. The value A^y g i should 
thus be subtracted from A7 in order to yield the correct in- 
trinsic shear. While this calculation gives an estimate of the 
shear signal due to contamination by lensed pairs, it does 
not take into account the (currently unknown) photometric 
redshift error distribution. In particular, E~ 1 [|z s — z;| < e] is 
close to zero because the lens-source separation is small, and 
hence its fractional systematic error may be large (possibly 
of order unity). 

The central values of the confidence intervals in Ta- 
ble 0] labeled "95% stat." are computed by subtracting the 
estimated gravitational lensing contamination from the ob- 
served intrinsic alignment signal (A7 — A7 9 ;). The statistical 
error on A7 9 ; is negligible and has not been included. 

The "99.9% stat.+sys." confidence interval includes, 
in addition to statistical errors, three of the possible sys- 
tematic errors: calibration biases in the shear measurement 
(Sec. 13.21 . errors in the "correction factor" R P /(R P — 1) 
of Eq. 1371 caused by seeing-induced correlations between 
high- and low-redshift galaxies (see Sec. I3.3.2L and sys- 
tematic errors in A7 9 ; (the statistical errors are negligi- 
ble) . The former is taken into account by dividing the com- 
puted shear by a worst-case calibration factor of 0.82 (ap- 
propriate if 57/7 = —0.18). The latter is taken into ac- 
count by assigning to R p a possible systematic error equal 
to 2|J? p [out] — 1|, where 7? p [out] is the pair ratio for the 
z a > zi + e sample. This is appropriate since ii p [out] = 1 
in the absence of magnification bias, photometric redshift 
errors, and seeing-induced systematics; it may even be a 
conservative estimate of seeing-induced systematics because 
some of the Tip [out] > 1 signal observed for brighter sub- 
samples is likely due to photometric redshift errors. Since 
the weakest constraints are always obtained from the small- 
est R p , we replace R p in Eq. I37H with R p — 2|i? p [out] — 1|. 
(In the outermost radial bin for the M r > — 18 luminosity 
subsample, this "minimum value" of R p is less than one, i.e. 
when systematic errors are taken into account we do not 
have a firm detection of clustering of source galaxies around 
this sample of lenses.) Finally, we allow A7 9 ( to range from 
zero to twice the measured value, i.e. the "maximum" value 
of the confidence interval for A7 is computed assuming no 
gravitational lensing contamination at all, and the "mini- 
mum" value is computed by subtracting 2A7 9 f instead of 
only A7 S ; from A7. 

As can be seen from Table 0| we find no evidence for 
intrinsic alignment of the satellite galaxies (there is a "detec- 
tion" at 95 per cent confidence for the M r < —22 subsample 
at 45 < r < IOO/1 -1 kpc, but given the large number of sub- 
samples investigated this should not be taken as evidence 
for A7 7^ 0). With all samples averaged together, we find 
conservatively that —0.0062 < A7 < +0.0066, although we 
cannot exclude the possibility of a large intrinsic alignment 
if present around only the faintest "lens" galaxies. 



Tangential Shear Comparison, Random Foregrounds 



- Z .i>urco >Z lo:i. ' 1 

Zi...-0- Kz <z,_..+ 0- 1 




Figure 13. The tangential shear signal computed around one of 
the random lens catalogs. 



5.2 Systematic shear 

Table includes limits on the systematic shear, i.e. the ad- 
ditive bias contributed by observational systematics to the 
total shear. Such a contribution could arise if, e.g. there were 
a correlation between the PSF and the survey boundaries. 
These limits were computed by taking the average value of 
the shears computed from the 78 random catalogs, and find- 
ing confidence intervals around this average using the Stu- 
dent's redistribution with 77 degrees of freedom (dividing a 
by \/78 because the measurement was done with 78 random 
catalogs). As shown, the systematic shear is consistent with 
zero for all samples; furthermore it is negligible compared 
to the statistical error in A7, as can be seen by comparing 
the upper limits in Table|S|to the statistical errors shown in 
Table 01 A plot of the systematic shear is shown in Fig. 1131 



6 DISCUSSION 

In this paper, we have set stringent new limits on intrinsic 
alignments as a contaminant of the galaxy-galaxy lensing 
signal. For example, averaging together all lens galaxies at 
radii out to 446ft" 1 kpc we derive -0.0062 < A7 < +0.0066 
(99.9 per cent confidence including identified systematics; 
see Table EJ or -0.0018 < A7 < +0.0021 (95 per cent con- 
fidence, statistical errors only). The statistical constraint is 
a factor of ~ 5 (~ 1.5 for our more conservative error bar) 
tighter than I A7I ^ 0.01 (95 per cent confidence) as o btained 
from APM/2dF data bv lBernstein fc Norberd" J2002I) . It also 
extends t o smaller radii th an the measurements of nearby 
spirals bv lLee fc Pen! <|200 lT) . Furthermore, our study of the 
luminosity and radius dependence of the intrinsic shear is 
useful for weak lensing studies that divide lens galaxies into 
luminosity subsamples. 

Our limits have been obtained by using "lens" galaxies 
with spectroscopic redshifts and "source" galaxies with pho- 
tometric redshifts, a different strategy than that selected by 
previous observers (iLee fc Penll200ll : iBernstein fc Norberd 
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Table 5. Limits on systematic shear in the source catalog 10 5 A7C; 
for various subsamples and ranges of radii. Values are quoted us- 
ing the errors computed from random catalogs (rand), at 95 and 
99.9 per cent confidence level, using the Student's t-distribution 
(with errors divided by v78 because there are 78 measurements 
of systematic shear). 



Subsample 


r 


10 5 A7< rami) 


M r 


h^ 1 kpc 


95% 


99.9% 


e = 0.1 


all 


[30,446] 


[-3.1,2.7] 


[-5.2,4.8] 


all 


[30, 100] 


[-4.6,20.2] 


[-13.5,29.1] 


aU 


[100, 300] 


[-5.6,1.8] 


[-8.2,4.4] 


all 


[300, 446] 


[-3.8,4.4] 


[-6.8,7.4] 


> -18 


[30,446] 


[-1.5, 13.1] 


[-6.7, 18.3] 


> -18 


[30, 100] 


[-13.3,59.9] 


[-39.5,86.1] 


> -18 


[100, 300] 


[-5.2,21.2] 


[-14.6,30.6] 


> -18 


[300, 446] 


[-6.2, 11.4] 


[-12.6, 17.8] 


-19 to -18 


[27, 493] 


[-6.9,3.9] 


[-10.8,7.8] 


-19 to -18 


[27, 100] 


[-14.3,40.9] 


[-34.1,60.7] 


-19 to -18 


[100, 300] 


[-9.5, 11.3] 


[-16.9, 18.7] 


-19 to -18 


[300, 493] 


[-11.0,3.4] 


[-16.1,8.5] 


-20 to -19 


[33, 545] 


[-5.8,3.0] 


[-9.0,6.2] 


-20 to -19 


[33, 100] 


[-21.7,25.3] 


[-38.5,42.1] 


-20 to -19 


[100, 300] 


[-10.0,5.0] 


[-15.3,10.3] 


-20 to -19 


[300, 545] 


[-6.1,3.9] 


[-9.7,7.5] 


-21 to -20 


[27, 735] 


[-4.6,2.4] 


[-7.1,4.9] 


-21 to -20 


[27, 100] 


[-17.7,29.7] 


[-34.7,46.7] 


-21 to -20 


[100, 300] 


[-11.9,3.9] 


[-17.5,9.5] 


-21 to -20 


[300, 735] 


[-4.3,2.9] 


[-6.9,5.5] 


-22 to -21 


[33,992] 


[-2.4,3.6] 


[-4.6,5.8] 


-22 to -21 


[33, 100] 


[-32.3,34.7] 


[-56.4,58.8] 


-22 to -21 


[100, 300] 


[-20.3,0.7] 


[-27.8,8.2] 


-22 to -21 


[300, 992] 


[-1.6,5.0] 


[-4.1,7.5] 


< -22 


[45, 992] 


[-9.7,7.7] 


[-16.0, 14.0] 


< -22 


[45, 100] 


[-120.7,72.7] [- 


-190.1, 142.1] 


< -22 


[100, 300] 


[-32.7,26.3] 


[-53.9,47.5] 


< -22 


[300, 992] 


[-8.6,8.2] 


[-14.7, 14.3] 


e = 0.05 


all 


[30,446] 


[-4.1,2.7] 


[-6.6,5.2] 


all 


[30, 100] 


[-5.4,28.2] 


[-17.4,40.2] 


aU 


[100, 300] 


[-5.9,3.9] 


[-9.4,7.4] 


all 


[300, 446] 


[-6.6,3.2] 


[-10.1,6.7] 


> -19 


[30, 446] 


[-3.7,8.1] 


[-7.9,12.3] 


> -19 


[30, 100] 


[-11.6,45.0] 


[-31.9,65.3] 


> -19 


[100, 300] 


[-2.9, 14.7] 


[-9.2,21.0] 


> -19 


[300, 446] 


[-9.9,6.3] 


[-15.6, 12.0] 


21 to -19 


[25, 545] 


[-7.8,0.8] 


[-10.9,3.9] 


21 to -19 


[25, 100] 


[-1.0,41.8] 


[-16.4,57.2] 


21 to -19 


[100, 300] 


[-11.6,2.0] 


[-16.6,7.0] 


21 to -19 


[300, 545] 


[-9.7,1.1] 


[-13.6,5.0] 


< -21 


[33, 992] 


[-5.1,4.1] 


[-8.5,7.5] 


< -21 


[33, 100] 


[-74.1,24.1] 


-109.3,59.3] 


< -21 


[100, 300] 


[-18.6, 10.0] 


[-28.9,20.3] 


< -21 


[300, 992] 


[-4.6,5.4] 


[-8.2,9.0] 



2002) who used lens-source pairs selected entirely on the ba- 
sis of spectroscopic redshift information. This strategy has 
dramatically increased our sample size and reduced our sta- 
tistical uncertainties for the brighter lenses (M r < —20) 
because these are typically at redshift z > 0.1, where most 
of their satellites are too faint to be targeted by SDSS spec- 



troscopy. This method does however have disadvantages. 
For our fainter-luminosity lenses, the photometric-redshift 
source method gives very large error bars (Table be- 
cause there are few physically associated lens-source pairs 
and hence any intrinsic alignment signal will be drowned 
out by the large number of unassociated pairs. The use of 
photometrically selected pairs also comes with the caveat 
that the limits on A7 could be circumvented by a small 
class of galaxies that have strong intrinsic alignments but 
always produce incorrect photometric redshifts. It is thus 
recommended that the values in Table 2] with redshift slice 
half- width e = 0.1 be used to constrain intrinsic alignments 
in galaxy-galaxy lensing studies, since in this case only about 
20 per cent of the physically associated source galaxies are 
scattered out of the sample by photometric redshift errors. 



Despite the recent interest in intrinsic alignments as 
a contaminant of cosmic shear (i.e. shear-shear correla- 
tions), there has been comparati vely little theoreti cal work 
on the density-shear correlation. iLee fc Pen! i200lT) do pro- 
vide a prediction for the density-shape correlation for spi- 
ral "source" galaxies, which we convert (see Appendix [Til 
to A7 = -0.004 at small radii r < l/i" 1 Mpc. This 
is marginally consistent with our results (it is within our 
more conservative confidence interval of —0.0062 < A7 < 
+0.0066), although the comparison is imperfect as the model 
applies only to spiral sources. More generally, measurements 
of the density-shear correlation A7 provide a constraint on 
intrinsic alignment models in addition to that provided by 
the shear autocorrelation measurements JPen et all l200fjl : 
iBrown et aT1l2002l) . It would thus be useful to compute 
the predictions of other intrinsic alignment models for A7 
and determine which models agree w ith the intrinsic shear 
measu rements presented here and in iBernstein fc Norberd 
J2002I) . 



We conclude by roughly estimating the intrinsic align- 
ment contamination to the gravitational lensing shear, (1 — 
J?~ 1 )A7. For example, if our source catalog restricted to 
z s > zi + 0.1 were used in a lensing study of the halos of 
SDSS spectroscopic galaxies, we would have a contamination 
fraction 1 - R~ 1 ~ 0.01 in the 100 < r < 300/T 1 kpc bin 
and intrinsic alignment contamination IA7I < 0.0095 (see 
Table 2J, corresponding to a contamination to the gravita- 
tional lensing shear at the level of < 9.5 x 10 . This is 
approximately equal to the size of the la error bars at these 
radii for our full sample (1.0 x 10 -4 ; see Fig, 111! and a factor 
of six less than the shear signal, ~ 5.8 x 10~ 4 . Slightly weaker 
constraints are placed on the intrinsic alignments around the 
brighter lenses, but in these cases the lensing signal is a fac- 
tor of several higher, e.g. for the —21 > M r > —22 magni- 
tude range we have | A-y] < 0.0124 at 100 < r < 300/i _1 kpc, 
1 — Rp 1 ~ 0.02, and shear signal jt = 2 x 10 -3 , so the con- 
tamination is at worst about a factor of eight smaller than 
the signal. These results indicate at 99.9 per cent confidence 
that, except possibly for the faintest lens galaxies, intrin- 
sic alignments are only contaminating the sub-Mpc scale 
galaxy-galaxy lensing signal at the < 15 per cent level. 
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APPENDIX A: INTRINSIC SHEAR AND 
CORRELATION FUNCTION 

In principle, the tangential component of intrinsic shear A-y t 
is a function not only of the transverse separation r between 
the lens and source galaxy, but also on the line-of-sight sep- 
aration u. The contribution of the intrinsic alignments to 
the shear -yt observed in galaxy-galaxy weak lensing is: 



7t(intr, r) = 



J Ajt (r, u)n(r, u)du 
J n(r, u)du 



(Al) 



where n(r, u)dud 2 r is the number of sources with line-of- 
sight separation between u and u + du in area d 2 r in the lens 
plane. This contamination is typically estimated by compar- 
ing the total number of sources observed per unit area in 
the lens plane, i.e. J n(r, u)du, to the number of sources ob- 
served in a randomly selected field, J n(oo, u)du. We express 
this in terms of a "pair ratio" R P (r): 



R P (r) = 



J n(r, u)du 
J n(oo, u)du 



(A2) 



The simplest interpretation of the pair ratio is that the frac- 
tion of the source galaxies that are physically unassociated 
with the lens is R ~ 1 , and that the remaining fraction 1 — R~ 
are physically associated with the lens. This motivates the 
definition of the projected intrinsic shear A-y(r) of physically 
associated galaxies as: 



Aj(r) 



7t(intr, r) 
1 — R v 1 



(A3) 



Most analyses of the effect of in trinsic alignments on galaxy- 
galaxy lensing stu dies (e.g. iBernstein fc Norberd 120021 : 
ISheldon et al]l2004T) have used Eq. 1A3I to compute 7t(intr) 
from the measured pair ratio R p and an estimate of the in- 
trinsic shear of physically associated galaxies A*y(r) based 
on spectroscopically selected lens-source pairs at the same 
redshift. 

While Eq. 1A3I is directly useful for measuring the in- 
trinsic alignment and estimating the contamination of weak 
lensing results, we can also write Ay(r) in terms of a projec- 
tion integral along the line of sight. Substituting Eqs. 1A11 
and 1A2II into Eq. IA3I . we find: 



Aj(r) = 



J A^t (r, u)n(r, u)du 
J[n(r, u) — n(oo, u)]du 



(A4) 



The integrands in both the numerator and denominator of 
Eq. IIA4(l vanish at large \u\ (so long as the magnification bias 
due to gravitational lensing is small). At small |u| (i.e. much 
less than the typical distance to the lenses, so that variations 
in n(oo,u) can be ignored), n(r,u) is related to the galaxy- 
galaxy correlation function via n(r, it) oc 1 + £(V r 2 + u 2 ). 
Therefore we may write 

J A 7t (r, u)[l+£{y/r 2 +u 2 )]du 



A-y(r 



J £(Vr 2 + u 2 )du 



(A5) 



which is the relation of A7(r) to fundamental (as opposed 
to observed) quantities. 
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APPENDIX B: RELATION OF INTRINSIC 
SHEAR TO POSITION ANGLE STATISTICS 



iLee fc Penl feOOCh and iLee fc Penl |200J) measure density- 
shape correlations using the position angle statistic, 



W2D(r 3 D) = ( cos2 



(cos 20), 



(Bl) 



where <f> is the position angle of the source galaxy's minor 
axis relative to the lens (0 indicates radial alignment, 7r/2 
tangential alignment) and the three-dimensio nal pair sepa- 
ration is r-io = \/r 2 + u 2 . ILee fc Penl (|200ll) consider spi- 
ral galaxies and hence interpret the minor axis as the (pro- 
jected) rotation axis of the galaxy, but Eq. (IBH is what is 
directly measured. 

We consider here a crude model for the conversion 
between u)2D and the intrinsic shear A7. Formulas for 
the transformat i on of ellipticity under shear are given by 
iMiralda-Escudel <ll99lf) : using that tan 20 = e^/et (in a 
coordinate system aligned with the direction to the lens 
galaxy) we can find the transformation of the position angle 
under shear, 



97* 



1 sin 2cf>. 



(B2) 



Substitution into Eq. (IBH yields 



du2D 



= \ /(-e" 1 81x120) ± cos 2^ = ^e" 1 ). (B3) 

We define this quantity to be the position angle responsiv- 
ity lZu>; for our r < 19 galaxies (with negligible measure- 
ment noise in the elliptic i ty), ( e -1 ) = 2.75 and TZ^ = 1.38. 
The theory of lLee fc Penl (2001) suggests that at small radii 
[r 3D < 1.5k' 1 Mpc) we should have U2d(x) = — 0.0057 for 
spirals (see their Fig. 2) or j t = u>2d/Hw — —0.0041. Since 
the galaxy autocorrelation function £(r3r>) 3> 1 at these 
radii, we expect A7 = —0.0041 fcf Ea. lX5l . This conversion 
is only a rough estimate since the relation Eq. lB3t between 
U2D and 7+ was derived for a gravitational shear, and it 
need not be the correct relation for intrinsic shear. 

Using the conversion = 1.38, the data presented in 



Fig. 2 of ILee fc Per] (1200 ID correspond to A7 
0.0025 at r 3D = 500/i _1 kpc. 



-0.0037 ± 



APPENDIX C: NOISE-RECTIFICATION BIAS 

Noise in galaxy images has several effects for weak lensing. 
The best-known effect is the introduction of measurement 
noise in the ellipticity. It has also been pointed out that noise 
can also couple to the PSF anisotropy to produce spurious 
power in cosmic shear me asurements, even if the PSF cor- 
rectio n is otherwise exact (lKaiserll2000l ; iBernstein fc Jarvisl 
2002). This occurs because of a noise-rectification bias, i.e. 
the ellipticity is a nonlinear function of the image inten- 
sity so that an unbiased estimate of the intensity need not 
translate into an unbiased estimate of the ellipticity. Here we 
estimate the effect of noise-rectification biases on the shear 
calibration of low signal-to-noise ratio galaxies. 

Our plan for approaching this problem is as follows: 

(i) Determine, for a circular Gaussian galaxy with unit 
radius a = 1 and unit central intensity 1(0) = 1, the mean 



and covariance of the measured adaptive moments m' to 
leading order in the noise. 

(ii) Use scaling and symmetry principles to derive the 

mean and covariance of K/l' ' for a general elliptical Gaussian 
galaxy. 

(iii) Propagate these errors to the PSF-corrected elliptic- 
ity e (/) . 

Here we will neglect corrections due to non-Gaussianity of 
the galaxy or PSF. 

We first introduce some notation. We write the nor- 
malization, centroid and moments of the galaxy as a six- 
dimensional vector: 



c a = (y4,x ,y , [M 



[M- L ] xy ,[M- L ] yy ) 



(CI) 



(We could have used the matrix elements of M to parame- 
terize the family of Gaussians instead of M -1 ; however using 
M _1 simplifies the algebra.) We now introduce the notation: 

r, a (r-r ) T M- 1 (r-r ) 
J(r;c )=Aexp . 



Thus the energy functional Eq. Q becomes 
1 



E{c a ;I) 



{I(r)-J(r,c a )Yd' 



(C2) 



(C3) 



or, written as an inner product: 

£( C V) = 1(/-J,I- J)- 

(We suppress the c a argument of J for clarity.) Minimizing 
this yields: 



(C4) 



= (I - J, J a ) . 



(C5) 



where J a = dJ/dc a . 

We may obtain the dependence of c a on I by taking the 
functional derivative of Eq. 1C5II with respect to I(ri): 



= J a (ri) + «/ - J, J aP ) - {J/3,J a )) 



5I(ri 



(C6) 



The second derivative of Eq. (IC5t is more complicated; for 
the particular case / = J, it is: 

dV 







SI(n)SI(r 2 

dc 13 dc' 



({J al 3,I-J)~{Jp,J a )) 
{{JaP, Jf) + (Jp-y, Ja) 



W(n) W(r 2 ) 

+ {J~/a, Jg) — {Jap-/, I — •/)) 

dc" r , s dc' 3 _ . , 
dl(n) ol(r 2 ) 

Defining H to be the matrix inverse Of [H" 1 ]^ = {Ja,Jp), 
and evaluating the derivatives at I = J, gives 



dc" 



SI(r) 
and 



= H a ?J a (r) 



2„/3 



d z c 



SI(r)SI(r) 



H ap ir* (J a ,J lS ). 



(C8) 



(C9) 



In Eq. (IC9II we have only computed the trace of the second 
derivative matrix; this is all we will need. 

A circular Gaussian galaxy with unit radius and central 
intensity has measured profile: 
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J(r)=. 



+ V (r) = J(r;C a )+v(r), 



(CIO) 



where 77 is the noise and C a = (1,0,0, 1,0, 1) are the mo- 
ments of the galaxy (in the absence of noise). The noise is 
assumed to be white: 

(V(ri)v(r 2 )} = N8 {2 \ ri - r 2 ). (Cll) 

The mean and covariance of the measured moments c a can 
then be determined to order 0{N) from the equations: 



(c a ) = C Q + I/V 



d c 



and: 



Cw{c a ,c?) 



■ N 



SI(r)SI(r) 

dc a dc" 



I(r) J(r) 



,2 

d r. 



Substituting Eqs. 1C8I and 1C9I yields: 



and 



-NH ap H 



rr^S 



Jf3{r)J lS {r)A r 



a0 



(C12) 



(C13) 



(C14) 



(C15) 



H=2 

7T 



(C16) 



The substitution of the particular functional form 
Eq. EH into Eqs. iCTH and (IHTst is a straightforward 
but lengthy exercise. The H matrix in the basis of Eq. lICH 
is 

( 1 1 1 \ 
1 

1 

1 4 
2 

\ 1 4 J 

and the product H lS J 1 s{r) appearing in Eq. 1C141 is 

H-> s J l6 (r) = -(p 4 - 2)e- p2/2 . (C17) 

From this we find that the mean values of the entries in M 1 
are 



<[M" 



c ) = {[M~ ] yy ) 



1 + V 

7T 



(C18) 



and {[M ^] xy ) — 0. Using the Taylor series for the moments: 



M 



(M 



(M 



0(M 



we can find the mean of M to be 

{Mii)= (i + ^-V) ft*. 

and the covariance to be: 

4 

Cov(My,M fei ) = -NiSikSjt + SuSjk). 



I) 3 , (C19) 
(C20) 

(C21) 



We now generalize Eqs. ljC2O0 and ljC210 . Clearly the 
energy functional is translation- invariant, so Eqs. )IC20IC21fl 
cannot depend on ro. If the amplitude A is increased by 
some factor, and the noise amplitude VJV is similarly in- 
creased, the energy minimization remains unchanged, so 
the result can depend on A and N only through the ratio 
Ajs/M . The problem is also invariant under linear transfor- 
mations of r, so that if M is changed, and M is increased 
by the factor %/det M (to account for the delta function 



in Eq. IC11I and the change of measure in the integration 
over d 2 r), then the linear transformation acts simply on 
Eqs. I|C20IC21(I . We have concluded that the generalizations 
of Eqs. HC20IC21H arc 



(Mij) = ( 1 + 4 
and 



M )M i3 
Tr^VdetM ' 



(C22) 



Cov(Mij,M k i) = 4 



AT 



7rA 2 VdetM 



(MikMji + MuMjk). (C23) 



The quantity Af /irA 2 ^ 'det M is recognized as v~ 2 , where v 
is the signal-to-noise ratio for det ection of the galaxy in a n 
adaptive elliptical Gaussian filter ^Bernstein fc ,Tarvisl2002^ . 

We have reached the final stage in our plan, namely to 
translate the errors in M into the calibration bias of the el- 
lipticity estimator. We consider here only the simplest case, 
that of a circular Gaussian PSF and an elliptical Gaussian 
galaxy with its major axis aligned along the x-axis. The es- 
timated ellipticity is: 



e (/) - 



)</) 



T(i) - T( p ) 



(C24) 



where T is the trace and Q = M xx — M yy is the + 
quadrupole moment. The effects of the bias (Eq. IC221 and 
noise (Eq. IC23I in M' 1 ' can be determined from the first 
and second derivatives (respectively) of efJ) with respect to 
M' 1 ^; the order 0(Af) terms are: 



(/), (SQ (1) ) Q {I) (5f^) 

e 1 -t- 



T(i) - r( p ) (tw - r( p )) 2 

Cov(QW,fW) Q (/) ((<5T«) 2 ) 
(T(i)-T( p )) 2 + {TO - T( p )) 3 ' 



(C25) 



where we have used the bias (5Q) = (Q) — Q. For a 
galaxy with only + ellipticity (i.e. position angle along 
the x- or y-axis), the biases and covariances obtained from 
Eqs. l|C22IC23(l are: 



{5Q (I) ) = 4u- 2 T 



<(<5T (/) ) 2 > = 4^- 2 T (I)2 (l + e« 2 ) 



Cov(Q (/) ,f (I) ) 
i(SQ (I) ) 2 ) 



= 8v~ 2 T W2 ef. 

= 4^- 2 r (/)2 (l + eV /)2 ) 



(C26) 



This leads us to the result for the change in ellipticity cali- 
bration: 



,(/) 



= —(l-3R^ + R- 2 + ey 



(/)2n 



(C27) 



We have verified Eqs. IH221 . JH2H . and llU27t in "toy" sim- 
ulations. 

We have used this result as an approximation to the 
shear calibration £7/7, although this is not exact because 
the shear calibration is determined by the derivative of 
{5e]P) with respect to ejvf , which is not exactly equal to 
the ratio {5e+)/e^P when the latter is not constant. The 
calibration bias obtained from Eq. (IC27II is shown in Fig. |S| 
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APPENDIX D: CHI-SQUARED TEST WITH 
ERRORS FROM SIMULATIONS 

We perform several tests for consistency with zero signal 
using a x 2 test. Since we have derived our covariance matrix 
from simulations, we must also take into account the noise 
from only having a finite number of simulations. 

If we have estimated a shear in N radial bins, we can 
construct an TV- dimensional shear vector 7. The M simula- 
tions provide us with M simulated shear vectors {"f^}a=i 
from which we can compute the sample covariance C (see 
Eg . 1331 . We are interested in the distribution of the variable 

X 2 =7 T r\ (Dl) 

under the null hypothesis that 7 is Gaussian with mean 
and the same covariance C as the simulated shear vectors. 

In order to determine the distribution of x 2 Monte 
Carlo, we define the vectors x = C _1//2 7 and x^ a ' = 
C — V 2 -^( a \ (Here the symmetric matrix C -1 ^ 2 has the same 
eigenvectors as C but the eigenvalues Xi are replaced by 
A~ 1/2 .) Then we find: 

X 2 = x T D- 1 x, (D2) 
where the matrix D is given by: 



D 



M-\ 



M 



E(a) <a)T - -1 
x x ' — xx 



(D3) 



Here x are the sample means of the as' . Since the Xi and 
xf* are independent Gaussians of zero mean and unit vari- 
ance, it is computationally easy to generate sample from 
Eq. I)D30 . so we can use a Monte Carlo method to compute 
the probability distribution for x 2 . 

As M —> 00 (with fixed N), D — * 1 and the distribution 
of x 2 converges to the x 2 distribution. However, the con- 
vergence is not rapid. To see this, we note that by Wick's 
theorem, D has expectation value equal to the identity 1 
and covariance: 

<(£>y - SvXDu - M> = SikS »+_ S fi k . (D4) 

Using the Taylor series for the inverse of a matrix, and the 
fact that the higher-order moments of the D distribution 
decrease as M~ 2 or faster, we find the asymptotic expansion: 

(D -1 ) = 1 + ((D — l) 2 ) + 0(M~ 2 ). (D5) 

Since x is independent of D: 

(X)=N+ N{N A + 1] +0(M- 2 ), (D6) 

which is found to agree with the mean obtained via Monte 
Carlo. Since the standard deviation of the x 2 distribution 
is V2N, we must take M > N 3/2 in order for the 0(1/M) 
correction term in Eq. 1D6I to be small compared to the 
width of the x 2 distribution. Therefore the standard x 2 test 
is only valid if M > TV 3 / 2 (not true for our random lens 
test, with M — 78 and N = 39). 



